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PREFACE 


The NSEG prograin was originally constructed by Mr. L. H. Leet of the 
United States Air Force Aeronautical Systems Division, Wright>Patterson 
Air Force Base. The code was subsequently modified and extended by 
Aerophysics Research Corporation under contract F33615-73-C-3039. The 
current version of the NSEG (Version III) extends the program appli- 
cability to higher speed (hypersonic turbo-ramjet) aircraft. It also 
includes various improvements generated by Mr. David T. Johnson of the 
Air Force Flight Dynamics Laboratory. The authors wish to extend their 
thanks to Mr. Walter Vahl of NASA for his extensive assistance during 
formulation and checkout of the turbo-ramjet propulsion system model now 
available in NSEG. The analytic basis of the turbo-ramjet model is 
due to Mr. Vahl. 

Mr. D. S. Hague of Aerophysics Research Corporation served as project 
leader for the present study. Dr. H. L. Rozendaal provided specialist 
support in the fields of propulsion system analysis and computer sciences. 
Mr. R. T. Jones, formerly of Aerophysics Research Corporation, has also 
made significant contributions to the NSEG code in studies preceding the 
present one. 

Additional details and copies of the program deck can be obtained from 
NASA Langley Research Center. 



NSEG - A SEGMENTED MISSION ANALYSIS PROGRAM FOR 
LOW AND HIGH SPEED AIRCRAFT 


D. S. Hague and H. L. Rozendaal 
AEROPHYSICS RESEARCH CORPORATION 

1 . INTRODUCTION 


Program NSEG is a rapid mission analysis code based on the use of approximate 
flight path equations of motion. Equation form varies with the segment type, 
for example, accelerations, climbs, cruises, descents, and decelerations. 
Realistic and detailed vehicle characteristics are specified in tabular form 
and a variety of layered atmosphere options are available. The mission 
specification is open-ended in that the upper limit on the number of flight 
segments to be included in a mission profile (currently one hundred and 
forty-nine) can be increased by increasing the size of a single common block 
(C0MTAB) above its current size of 3000 words. The code contains an English 
language oriented input procedure for describing the mission segment sequence 
to be employed. In addition to its mission performance calculation capa- 
bilities the code also contains extensive flight envelope performance mapping 
capabilities. For example, rate-of-climb, turn rates, and energy maneuver- 
ability parameter values may be mapped in the Mach-altitude plane. Where 
suitable graphics capabilities exist these maps may be drawn by machine in 
the form of contour plots. 

The code contains several approximate flight path optimization capabilities 
based on Rutowski energy-like criteria. These flight path optimization 
formulations permit inclusion of minimum time or fuel flight segments 
and maximum range segments during climb or descent segments. Approximate 
take off and landing analyses are also performed. At high speeds centri- 
fugal lift effects are accounted for. Extensive turbojet and ramjet engine 
scaling procedures are incorporated in the code. Take off and landing 
analyses are also available which employ a high lift aerodynamic analysis 
model based on the Air Force Flight Dynamics Laboratory DATCOM method. 
Alternatively, user supplied high lift aerodynamics can be employed. 

This report is Volume II of 3 volumes. Total program documentation consists 
of: 

Volume I. Theoretical Development 

Volume II. Program User's Manual 

Volume III. Test Problems 



2 . PROGRAM STRUCTURE 


NSEG has eight overlays. Program structiore is defined in Figure 1-. 

The main control program is NSEG II. This program calls the next level of 
overlay, 1.0, 2.0, and 3.0 directly. Overlays 3.1 and 3.2 are, in turn, 
directly called from the program MISION of overlay 3.0. The remaining two 
overlays, 3.3 and 3.4, are called from subroutine FLITEN which is also 
directly called from program MISION. 

Gross functions of each overlay are as follows: 

OVERLAY 0.0, NSEG II - Main controlling program. 


OVERLAY 1.0, INISEG - Initialize the program including common 
blocks and input directories. 

OVERLAY 2.0, INSEG - Read the input data. 

OVERLAY 5.0, MISION Control program for mission and performance 
calculations . 


OVERLAY 3.1, QUTATM - Controls printout of the atmospheric properties 
employed. 


OVERLAY 3.2, ENSEG - Controls generalized performance calculations. 
OVERLAY 5.3, FLIRAN - Generates data for performance contour maps . 
OVERLAY 3.4, FLIMAP - Produces performance contour maps. 


The program consists of approximately ten boxes of cards (twenty thousand source 
cards). Coding is predominantly written in the algebraic FORTRAN language. 

The program currently operates on the CDC 6400 and 6600 series computers. 

Core requirement depends strongly on the computer center at which the 
program is run. At Lawrence Radiation Center the retractable binary loads 
are IISOOO^, at the Control Data Cybernet Centers 135000 q, and at Langley 
Research Center 150000_. 

O 

Constants utilized in the NSEG code employ English Units (lb., ft., sec., naut. 
mi., and degrees Rankine) exclusively. It is imperative therefore that input 
quantities employ English Units; constants utilized in the illustrations 
reproduced in the text reflect this convention. Provisions of NASA Policy 
Directive (NPD 220.4) have been waived for those portions of this report that 
pertain to the NSEG computer code. 
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3. SUBROUTINE LIST 


AHNACH 


ALTA 


ASCALE 


ATNOSD 


ATM053 


ATM0S9 


ASROi 


AUGHT 


Writes out the z array which contains scaled engine data vs. 
Mach number and altitude. 


Entry point in ATMOS D 


Constructs the z array from the atmosphere routine 


Controls atmospheric property calculations 


1962 standard atmosphere (Original NSEG) 


Arbitrary ten 5 >erature profile atmosphere 


Looks for a table location associated with a given starting 
location in common block /TABNAM/ which contains names and 
locations of input table names and data, respectively. 


Dusony routine to supply thrust augmentation 


BCDDEC Part of input procedure 



BCDINT 


BERAN 


BSCALE 


CBFUEL 


CHANGE 


CHECK 


CHMACH 


CHRNUM 


CLIMB 
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Part of input procjedure 


Calculates altitude. for maxiimiiB range factor given Mach number 
and weight 


Assists in setting up the Z array (see AHMACH) from the atmosphere 
routine 


Computes time to turn through specified angle; in Option 17 
Combat 


Sets NSR to zero 


Provides instantaneous conditions when lift equals weight 


Writes out the z array (see AHMACH) 


Detects illegal input characters on data cards 


Calculates climbs and descents 



COMP AC Compacts card images 


CMAX Calculates Mach number for maximum rate of climb or minimum 
fuel climb given altitude and weight 

CONPLT Contour plotting routine for performance maps in overla>^ 3.4 

CRCOMP Computes range factor optional 

CSCALE Sets up atmospheric z array 

DBOUT Sets up tape 25 

DISC0T Two dimensional table lookup subroutine 

OORDER Sorts input list in alphabetic and numeric order 

DRAG Aerodynamic control routine 

DSERCH Carries out directory search for input 
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DUMPRG 


DUMPTZ 


EDMEIPT 


EDCLMB 


EN/VX 


ELT 


ENCASE 


ENDAT 


ENERGY 


ENGINS 


Duaps the RG array > the array which contains data relative to 
completed segments (e.g. option, range, weight, fuel, etc.) 


Dumps the TZ array* the array which contains data necessary for 
the serial execution of flight path segments as specified by 
input data 


Locates maximum performance criteria along specific energy contour 


Performs maximum E cliiiA> 


Calculates Mach number for maximum endurance given weight and 
altitude 


Provides elaosed computer time computation 


Controls automatic selection of propulsion option. 


Calculates thrust and specific fuel consumption 


Computes energy maneuverability parameters at specified G-level, 
Mach number, weight, and altitude 


Engine simulation subroutine for the turbojet, ramjet and/or 
combined engine propulsion option. 
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ENSEG 


ContTote mtes-ion ‘eimutaPions 


(ETogram) 


ENSET 


ERROR 


FCLMAX 


FLYMHP 


FETCH 


FILLPM 


FI LBCD 


Modifies input propulsion data basis for tabular thrust and 
fuel flow option 


Reports table .errors to user 


Finds Cl Max 


Flies the vehicle between two M-h points 


Reads each input card 


Initializes vehicle data mission segment options 


Part of input routine 


FNLOSS 


User supplied routine to given propulsion losses 


FUELIF 


FILLTZ 


FLIMHC 


FLITEN 


FLniAP 


FLIRAN 


FLIENV 


FLIXEP 


FLIXPR 


Computes in-flight fuel used j its disposition « tank':drop tinfts, 
and fuel remaining 


Fills the C array from the term array and gives optional print 
of term array. The term array is a 10 word arra;y of segment - 
related data created by the input option selected (INXX1-INSS21) 
The C array contains the term arrays in serial fashion for 
subsequent mission execution. 


Constant altitude-constant Mach number cruise 


Controls generalized performance calculations 


Controla perforrnance plots 


(Program) 


Controls generalized performance calculations 


(Program) 


Given weight and a set of Mach numbers, computes flight envelope 


Computes endurance performance 


Begins range performance calculations 
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FLIXCP 


Computes climb perfonnance 


FLITOL 


FLIMLD 


FUELLD 


GTI1I2 


GOGO 


HDAMPT 


HDCLMB 


HIHO 


Controls take-off and landing calculation 


Cong)utes range factor performance 


Computes initial fuel load, its disposition, and useable fuel 
load 


Part of table look up routine 


Calculates constant Mach or C^ cruises where independent variable 
is time or distance 


Finds maximum of performance criteria at constant altitude 


Flies maximum rate of climb path 


Controls N dimensional table look up routine 
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HNSEG 


INIMHP 


INDXCK 


INSEG 


IVNAME 


INISEG 


INITZ 


INIWGT 


INXOl. . 
INX22 


LANDIN 


LANDNG 


Utility routine; prints out a heading identifying the NSEG 
program version 


Initialize M-h point calculation 


Detects incorrect data input names 


Controls error check in data input (Program) 


Decomposes data input cards 


Initializes USEG II (Program) 


Initializes the TZ array, the array which contains the data 
necessary for the serial execution of flight path segments as 
specified by input data 


Initial weight breakdown computations 


Input routines for Options 1 through 22. Written in Macro-Foftran 
language and assembled on GDC 6600 Fortran. The present program 
uses the GDC 6600 Fortran versions directly without the need to 
use Macro-Fortran 


Provides organized print out of landing inputs 


Performs landing calculations 
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LINES Controls the niunber of lines printed on each page 


LISTZL Notates a plot 


LNCU4B Controls linear Mach-altitude path calculation 


MAPOl Prepare E map 


MAP02 Prepares E/m map 


MAP03 Prepares L/D Map 


MAP04 


Prepares range factor map> 



MAP05 prepares thrust map 


MAP06 Prepares drag map 



MAP07 


Prepares specific fuel consumption map 


M\P08 Prepares fuel flow map 


MAP09 Prepares specific energy map 


MAPIO Prepares lift/ (thrust-drag) m^ 


MAPll Prepares radius of turn map 


MAPI 2 Prepares time to turn map 


MAPOUT Controls map outputs 


MDCLMB Flies best acceleration path 


MDHPT Finds maximum performance criteria at constant Mach number 
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MHPTAB 


Outputs Mach-altitude path histories in tabular form 


MHPATH 

MISION 

NDATAl 

NDATA2 

NDATA3 

NDATA4 

NDTLV 

OBSHOK 

OPTXXI 

OPTXX2 


Draws a Mach- altitude path 


Main program for Overlay 3, which performs all flight path 
and flight envelope calculations 


Directory of input data 


Directory of table names 


Extension of input data directory 


Contains data (temperature and pressure vs. altitude) to 
describe an arbitrary atmosphere used in ATM0S9 


N-dimensional table look-up routine 


Given supersonic flow properties ahead of a wedge this routine 
computes flow properties behind the resulting shock wave. Two 
dimensional isentropic flow is assumed 


Performs Option 1 calculation, climb or descent path 


Performs Option 2 calculation, maximum rate of climb or minimum 
sink rate 
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0PTXX3 


0PTXX4 


0PTXX5 


0PTXX6 


0PTXX7 


0PTXX8 


0PTXX9 


OPTXIO 


Performs Option 3 calculatiopf Radius mission option 


Performs Option 4 calculation,. cruise, climb tO; specified weight 


Performs Option 5 calculation, cruise climb for specified range 
or time 


Performs Option 6 calculation > constant altitude cruise to specified 
end weight 


Performs Option 7 calculation, constant altitude cruise for’ a 
specified range 


Performs Option 8 calculation, constant altitude cruise for a 
given time 


Performs Option 9 calculation, buddy refuel 


Performs Option 10 calculation, transfer Mach-altitude weight 
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OPTXll 


OPTX12 


OPTX13 


0PTX14 


OPTX15 


OPTX16 


OPTX17 


OPTX18 


Perforas Option 11 calculations, NSEG segment selection logic 


Perforns Option 12 calculation, instantaneous weight change 


Perforas Option 13 calculation, instantaneous Mach-altitude change 
at constant wei^t 


Perforas Option 14 calculations, restores previous state or 
computes point calculation 


Perforas Option 15 calculation, adjust range in previous segment 
to give specified total range 


Performs Option 16 calculation, cliab/accelerate between specified 
N-h points 


Perforas Option 17 calculations, fuel weight change in warm up, 
take-off, loiter, or coirt>at 


Perforas a store drop 
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0PTX19 


OPTX20 


0PTX21 


OPTX22 


OUTATM 


OUTXXX 


OVWGTE 


PAGHED 


PATH 


PAYLOD 


Performs a tank drop 


Performs Option 20 calculation, fuel allowance for specified 
time at given power setting or T/W ratio 


Performs engine scaling when dual engine propulsion is selected 


Range matching generalized iterative control option 


Atmosphere output control program 


(Program) 


Provides output from mission segment options 


Computes the overall weight empty including stores, tank, 
and pylon combinations 


Provides the output page heading 


Computes conditions along a climb or descent arc 


Computes vehicle payload including store, tank, and pylon 
contributions 
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PERPLT Control program for the twelve performance maps 


PHMACH Provides organized print of Z array (see AHMACH) 


PICINX Mission segment input control routine 


PLWDRP To be used for payload drop 


PNAMES Writes directory of input names 


QCLIMB Computes constant dynamic pressure climb 


QBOUND Imposes dynamic pressure boundary limits on plots 


RETRNW Computes aircraft return weight 


RFMAX Computes Mach number and altitude combination for best range 

factor given weight; or Mach number for best range factor given 
altitude and weight 

17 
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RGSAVE 


Saves the RG array (see DUMPRG) ; the RG array will be 
printed if IPRGSV is nonzero 


RNSEG Input main control program 


RPOLY Computes the roots of a polynomial, given the polynomial coefficients 


RUN Calculates cruise paths 


SACSP Store and pylon aerodynamics 


SACTP Tank and pylon aerodynamics 


SACNSl NSEG Option 1 clean vehicle aerodynamics 


SACNS2 NSEG Option 2 clean vehicle aerodynamics 


SEGINI Initialize segment calculation 


SEGOUT Outputs segment calculations on Unit 6 and Uiit 15 in obser- 
vation function style 
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SETYIC 


Initializes Y, Yl, Y2, Y3 arrays 


SHELL Sorts array o£ names into algebraic order 


SHELLX Sorts array of values into numeric order 


SLVALU Part of table look up package 


STASH Utility routine in input procedure 


STORES Computes the number of stores dropped and the dropped weight 


STORWT Computes total store drop and return weights 


TABRE Controls table data read 


TABUPl Part of table look up routine 


19 



TAKEOF 


Prints take-off data input 


TAKOFF Carries out take-off calculation and outputs result 


TANKS Computes current ntunber of tanks and pylons dropped and their 
drop weight 


TNKSWT Confutes total tank drop and return weight 


TLAERO Take-off and landing aerodynamics option 


TLUl Two-dimensional table look up 


TOP Calculates ceiling for a given rate of climb and weight 


TLU Two-dimensional table look up routine 
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TSRCH 


TZRGDB 


UNPACK 


UNIT25 


VMAX 


VOIDXX 


XALT 


XCEL 


XCLECT 


Compiles starting locations for data tables in common block 
C0MTAB and stores these locations in common block TABI^AM 


Writes TZ and RG arrays on Unit 25 


Unpacks the MTZ array and writes the results and the TZ array 
on Unit 6 


Controls input and output for generalized performance calcu- 
lations 


Calculates maximum Mach number for steady level flight given 
altitude and weight 


Sets weight, altitude, and Mach number from RG array 


Computes altitude for a given pressure ratio 


Calculates constant altitude accelerations and decelerations 


Selects the next mission segment option and passes on the TZ 
and RG arrays 



XFIT2 Second order curve fit routine 


XNLTST Searches card image for non-numeric character 


XTRACT Decoo^oses card images 


XXINIT Initializes the XX array 
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4. DATA DECK STRUCTURE 


The basic NSEG data deck is presented in Figure 2. Major input blocks are 

1. Three TITLE cards . The three TITLE cards will be used 

as a heading and will appear on all regular NSEG , 

print out pages along with a page number. 

2. STCASE/TABSIZ/ . This data set is used to reserve 
space in the internal table array for the various 

tables to be used during the computer run. ; ^ . 

3. STCASE/DATAl/ . This data set is used to speciJfy 
values for all vehicle characteristic tabular data 
and data used to control the print output of the 

program. c , , 

4. STCASE/GENPC/ . This data set is used to specify 
data for the general performance calculations and/or 
performance maps . This data set is optional and 
need only be included if general performance calcu- 
lations and/or performance maps are to be generated. 

5. S TCASE/MSEG/ . This data set is used to specify details 
of the mission to be flown. 

6. RUN . The run card is the indicator to the NSEG 
input package to stop reading data and execute the 
t£isk defined by the data just processed. 

Additional data cases may follow the run card; however, the following deck 
structure must be observed: 

1. Each case must start with three title cards. 

2. The STCASE/TABSIZ/ data set need not be repeated 
in subsequent cases unless the table sizes are 
to be altered. If the STCASE/TABSIZ/ data set 

is included in a subsequent case, then all tabular 
data must be input in that case since the starting 
locations for the various tables were redefined by 
the STCASE/TABSIZ/ data set. 

3. The STCASE/DATAl/ data set need not be repeated in 
subsequent cases unless data items in that set are 
to be altered. It should be noted that all vari- 
ables will have the same value as established by 
the previous case unless they are changed by a data 
card in the current case. Therefore, only data 
items whose values are to be changed need be given 
in the current case. 
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d. The STCASE/GENPC/ data set is destroyed at the end of 

each case; therefore, if general performance calculations 
are required in a subsequent case, all data necessary to 
define the desired performance calculations must be 
input each time. 

Some typical exaII^)les of additional NSEG data sets are presented in 
Figure 3. 
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5. ATMOSPHERE OPTIONS 


NSEG contains the original NSEG 1962 U. S. Standard Atmosphere and an 
arbitrary layered atmosphere in which the ten^)eratures or temperature and 
pressure at the base of each layer are input. Aerophysics Research 
Corporation has the following atmospheric models available in NSEG 
format on request. 


1. ATMOSl: 1962 U. S. Standard Atmosphere, from AFFDL 

Gentry Arbitrary Hypersonic Body Program 

2. ATMOS 2 ; 1962 U. S. Standard Atmosphere, from Langley 

Research Center STOP program for supersonic transport 
optimization 


3. 


4. 


5. 


6 . 


7. 


8 . 


9. 


ATM0S3: 1962 U. S. Standard Atmosphere; analytic 

model from original NSEG 

ATM0S4: 1963 Patrick Air Force Base atmosphere 

ATMS59 : 1959 ARDC atmosphere 

ATJAN6: 1966 January atmosphere at 30® North 

ATJUL6 : 1966 July atmosphere at 30® North 

ATMS62: 1962 atmosphere from AFFDL ATOP program 

ATM0S9: arbitrary atmosphere 


As noted above, only atmospheric models 3 and 9 are now an integral part of 
NSEG. The other atmospheric models can be added and called from ATMOSD 
by a simple progreun addition. 


IPATMD 

INDATM 


5 . 1 Data Input for Atmosphere 

= 0, Don’t print out atmosphere model 
= 1 , Print out atmosphere model 

= 3, Use ATM0S3 
= 9, Use ATMOS9 


When INDATM 


NHG 

HGi 

TM^ 


-9, the following additional data is required: 

= Number of layers; nominally not to exceed 25 
= Array of altitude values , in feet 
s Array of temperature values, ® F absolute 


2 

» Array of pressures, lb. /ft. . NOTE that pressures can be 
detemined from sea level pressure and temperature profile. 
In this case set 

HG^ = -16505. 

* Desired temperature profile 
PGi = 3711.0839 


26 



6. PROPULSION 


6.1 Simplified Propulsion Systems - Introduction 

The simplified propulsion options in NSEG provide for the computation 
of a vehicle's thrust force and fuel consumption primarily by table 
lookup procedures. Any one of three propulsion systems, dry, wet, or 
maximum, may be employed at any instant along the flight path. Provision 
is made for automatically switching from one propulsion system to another 
on Mach number. Alternatively, the user may directly specify which of 
the three propulsion systems is to be employed. 

In all options the thrust force, Fjij, is specified in the form 

= FN£ (h , M) ; i = 1 , 2 , 3 

In the first two options the fuel flow rate 

FF = FF(PS, M, h); i = 1, 2 

In the third option 

FFj = FFj(h, M) 

In the automatic Mach number propulsion system switching mode, the following 
sequence is employed 

F^^, FF^; M^l.O 

Fj^, FF = FN3, FF3; 1.0 < M ^ 3.5 
Ff^2» FF2; M > 3.5 

Within the program the propulsion system input is scaled by appropriate 
factors to produce the following modified thrust and fuel flow data 

(Fn)! = FNi/DR ; i = 1, 2, 3 

(FF)^ = FF^/(100 X TR^'^) ; i ^ 1, 2 

(FF)^ = FF3/TR^*^ 

where DR is the atmospheric density ratio and TR is the atmospheric temperature 
ratio. 

In the dry and wet options (i = 1 and 2, respectively), the power setting is 
specified as a percentage of maximum; hence, the factor of 100 in the scaling 
equations for fuel flow in the cases i = 1 and 2. Again, in the dry and wet 
options, three power setting options are possible. These correspond to 

Maximum available in system 
Fj^ * Required to equal drag 

Specified value of power setting 
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The corresponding fuel flows are determined from the power setting 


PS = FN/(FN Maximum Available) 

It should be noted that power setting data is specified in the mission 
segment options of Seations 9,1 - 9,20 and Iwnoe does not form part of 
the basic propulsion system input. 


6.1.1 Relevant Subroutines 

The subroutines of particular significance ir the simplified propulsion 

system calculations are listed below, 

AHMACH - Used to write out the input and scaled thrust data in an 

organized manner. It is called from the propulsion system 
scaling routine ENSET. 

ASCALE - Used to scale, the input thrust data. All input thrust data 
for the dry, wet, and maximum thrust options are divided by 
density ratio. This data resides in TTABOl, TTAB03, and 
TTAB05. Called from ENSET. 

BSCALE - Used to scale the input fuel flow data in the dry and wet 
options. All input fuel flow data in TTAB02 and TTAB04 
is scaled by dividing by (100 x temperature ratio ** 1.5). 
Called from ENSET. 

CSCALE - Used to scale the input fuel flow data in the maximum power 
option. Input fuel flow data in TTAB06 is divided by 
(temperature ratio ** 1.5). Called from ENSET. 

CHM^CH “ Used to write out the input and scaled fuel flow data for the 
maximum thrust case in an organized manner. Called from 
ENSET. 

ENCASE - Controls automatic propulsion system switching on the basis 
of Mach number. Called from ENDAT. 

ENDAT - Performs the actual thrust and fuel flow calculations. It 
is called from a variety of mission segment and point per- 
formance calculations. 

ENSET - Controls the initial propulsion system input scaling and 
optional print out. 

PHMACH - Used to write out the input and scaled fuel flow data for 
the dry and wet cases in an organized manner. Called from 
ENSET. 

Flow charts for ENCASE, ENDAT, and ENSET are presented in Figures 4, 5, 


and 6. 
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6.1.2 Data Input 


6 . 1 . 2 . 1 General 


SYMBOL 

DESCRIPTION 


CASEN 

Engine option indicator: 

= 1, dry 
- 2 , wet 
=3, maximum 

= 4, automatic selection on Mach number 


FFLMIN 

Minimum fuel flow 


IPENDR 

Raw engine data print indicator: 
= 0, omit raw data print 
= 1 , print raw data 


IPENDS 

Scaled engine data print indicator: 
=0, omit scaled data print 
= 1, print scaled data 


NUMENG 

Number of engines 


PSET 

Power setting indicator: 

= 1, maximum thrust in CASEN = 1, 2, or 
= 2 , minimum thrust in CASEN = 1 , 2 , or 
= 3 , specified power setting in CASEN = 
or 3 

3 

3 

1, 2, 

SFCFAC 

Multiplier for SFC value 


TSIZEF 

Engine scale factor 


THRMIN 

Minimum thrust 


6 . 1 . 2 .: 

Dry Power, CASEN = 1 


SYMBOL 

DESCRIPTIC»i 


INDTOl 

Thrust table indicator: 
= 0, no thrust 
= 1 , thrust data given 


TTABOl 

Maximum dry thrust: T = T(h, M) 


ITOIX 

Number of altitude values in TTABOl 
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SYMBOL 

DESCRIPTION 

ITOIY 

Number of Mach number values in TTABOl 

INDT02 

Fuel flow table indicator: 
= 0, no fuel flow 
= 1, fuel flow data given 

1TAB02 

Dry thrust fuel flow Table in form of S.F.C. 


W = W(PSD, M, h) 


where PSD is per cent maximum 

IT02X 

Number of PSD values in TTAB02 

IT02Y 

Number of Mach number values in TTAB02 

IT02W 

Nximber of altitudes in TTAB02 


6. 1 . 2 . 3 Wet Power, CASEN = 2 


SYMBOL 

DESCRIPTION 

INDT03 

/ 

Thrust table indicator: 

= 0, no thrust 

=1, thrust data given 

TTAB03 

Maximum wet thrust: 


T = T(h, M) 

IT03X 

Number of altitude values in TTAB03 

IT03Y 

Number of Mach number values in TTAB03 

INDT04 

Fuel flow table indicator; 
= 0, no fuel flow 
= 1, fuel flow data given 

TTAB04 

Wet thrust fuel flow table in form of S.F.C. 


W = W(PSD, M, h) 


where PSD is per cent meucimum 

IT04X 

Number of PSD values in TTAB04 

IT04Y 

Number of Mach number values in TTAB04 

IT04W 

Number of altitude values in TTAB04 
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6. 1.2.4 Maximum Power » CASEN = 3 


SYMBOL 

DESCRIPTION 

INDT05 

Thrust table indicator: 
= 0 , no thrust 
s 1, thrust data given 

TTAB05 

Maximum power thrust table 


T « T(h, M) 

IT05X 

Number of altitude values in TTAB05 

IT05Y 

Number of Mach number values in TTAB05 

INDT06 

Fuel flow table indicator: 
» 0, no fuel flow 
= 1, fuel flow data given 

TTAB06 

Maximum power fuel flow tablein form of S.F.C. 

W = W(h, M) 

(Max. PS data only) 

IT06X 

Number of altitude values in TTAB06 

IT06Y 

Number of Mach number values in TTAB06 


6. 1.2. 5 Automatic Propulsion Switchings CASEN » 4 


When CASEN « 4.0, the propulsion option is internally switched among CASEN = 1, 
2 , and 3 according to Mach number as follows; 

Mach CASEN 

0 - 1.0 
1.0 - 3.5 
> 3.5 


1.0 

3.0 

2.0 





6.2 Turbojet, Ramjet and/or Combined Engines 

with Inlet. Precompression Effects (Combined Engine Option) 


The turbojet, ramjet and/or combined engine simulation option provided in NSEG 
by the ENGINS subroutine provides a more realistic engine simulation 
option than those described in Section 6.1 due to the following features: 

1. Inlet flow field precompression effects due to airstream 
flow deflection are accounted for in thrust calculations. 

2. A more accurate description of engine performance is provided 
by tabular engine data. (see 6.2.3) 

3. Turbojet only, ramjet only, or turbojet /ram jet combination 
engine options are provided; in-flight option selection as a 
function of Mach number can be employed. 

4. Engine-related lift and drag forces are computed and accounted 
for in the net thrust and the desired vehicle lift. 

5. High altitude and Mach number atmospheric properties related 
to engine performance are described by tabular data. 

The following sections provide an analytical description of the turbojet, 
ramjet and/or combined engine propulsion option and complete user-related 
information. 


6.2.1 Analytical Description 

The basic assumption of two-dimensional flow is employed in the turbojet, 
ramjet and/or combined engine simulation subroutine. The geometry and 
various flow regions relevant to the simulation are shown in Figure 7. 

The computations performed are divided into seven main categories as follows: 

1. Free stream properties at Station 0 which are a function 
of aircraft Mach number and altitude. 

2. Flow field conditions at Station 1, aft of the wing shock. 

3. Inlet computations which yield the turbojet inlet recovery 
ratio and ram drag. 

4. Flow field conditions at Station 2, aft of the engine inlet 
wedge shock. 

5. Turbojet thrust, airflow, and fuel flow computations, a 
function of flow field conditions at Stations 1 and 3 (at 
the turbojet compressor face) and throttle setting. 

6. Ramjet thrust and fuel flow computations, a function of 
flow field conditions at Station 2 and M^, 
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7. Spillage drag computations, if any, are performed either 
in the turbojet iteration loop for the turbojet alone 
engine option or in the ramjet iteration loop for the ramjet 
alone engine option, or the combined engine option. 

Details of these computations are described in the following text. 

6. 2. 1.1 FREESTREAM PROPERTIES 


For freestream conditions where M ^ 3,5 and altitude exceeds 19812 meters 
(bSOOO ft} real gas freestream total pressure (pto) total temperature 
values are computed by the two-dimensional table lookup subroutine DISCOT. 
Linear interpolation is utilized. Plots typical of these data appear in 
Figure 9. Freestream total pressure in p.s.f. is then given by 


P = 2116. 22 *p* (p* in atmospheres) 

0 

If the above Mach-altitude conditions are not satisfied, total temperature 
and total pressure values are computed by first determining freestream static 
conditions (PqjTq) from the 1962 Standard Atmosphere Subroutine ATMS62. 

Static enthalpy (H) and pressure ratio (pro) then determined from 

Table II, given the static temperature. Total enthalpy (Hf-) can next be 
computed by the expression 

Hj, = H + (Mq * aQ)^/50073.2 

Given total enthalpy, the freestream total temperature and total pressure 
ratio (prt ) determined from Table II. Freestream total pressure 

is then ° given by 


'0 




6. 2. 1.2 FLOWFIELD CONDITIONS AT STATION 1 


Given the free-stream flowfield parameters at Station 0, their counterparts at 
Station 1 aft of the wing shock are determined by the two-dimensional 
isentropic shock relations (NASA Report 1135) . These relationships are 
mechanized in the subroutine OBSHOK which performs oblique shock compu- 
tations for wedge angles below the critical value and normal shock 
computations at higher wedge angle values. 

6. 2. 1.3 INLET AND RAM FORCE COMPUTATIONS 


Given flowfield parameters at Station 1 and the engine streamtube area, the 
ram drag is computed from the equation 

TURBOJET RAM FORCE = * A.J, * (M^ * a^)^ + A.J. (P^ - P^) 

where A,^ is the area of the stream tube of the airflow utilized by the 

turbojet. Ramjet ram forces are accounted for a priori in the ramjet 
specific impulse data. The ram force computed above, parallel to the 
vehicle wing under surface, is subsequently resolved into components in 
the lift and drag directions. 
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The engine capture ratio CA /ApuL^] is determined from Figure 9, 

Table III, given The air flow rate captured by the inlet then 

becomes 

^A^ "" ^ApuLL * ^'^c^^FULL^ " ^1 * ^INLET * ”1 * * CA^/ApyLL^ 

Capture ratio less than 1.0 indicates operation at a below design Mach 
number. As seen in Figure 7 the result of below design operation is that 
the wedge shock fails to intersect the inlet cowl resulting in inlet air 
spillage. This spillage produces forces (inlet spillage lift and drag) 
which are computed using a combination of continuity relationships and 
geometry as shown in Figure 8. 

6. 2. 1.4 FLOWFIELD CONDITIONS AT STATION 2 


Given the flowfield parameters at Station 1, conditions at Station 2 
are computed using subroutine OBSHOK as in Section 6. 2. 1.2 above. That 
is, the two-dimensional flow behind the precompression (wing) surface shock 
is turned again through the inlet wedge angle (6j). 

6. 2. 1.5 TURBOJET AIRFLOW, FUEL FLOW, AND THRUST COMPUTATIONS 


For a base size turbojet, airflow and fuel flow requirements at full 
throttle are first determined via the one-dimensional table lookup 
routine FTLUP. The corresponding maximum thrust is determined via the 
two-dimensional table lookup routine DISCOT. The corrected airflow (W^/6/6) 
requirement is determined from Table V, given the total temperature at 
Station 3, Tt?' This value is then used to determine the turbojet airflow 
at any flight “^condition by the equation 


where 


W, 


(W^) 

corrected 


_6 

/e 


* Ft- 


(W^) 

corrected 


the base engine corrected airflow rate 


6 


/2116.22» ratio of compressor face total pressure to 
3 sea level static pressure 


0 “ T^ /518.67, ratio of compressor face total temperature to 

3 sea level static temperature 

F-fj = the turbojet scaling factor obtained by specifying 

(1) desired net thrust at sea level standard conditions 
or (2) desired turbojet airflow at sea level standard 
conditions 
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Note; It is assumed that Tt, = T^ = Tt, = Tt^ (Figure 7). 

^210 



The fuel flow rate to air flow ratio is determined from Figure 9, Table VII, 
given The maximum turbojet thrust to air flow rate is found from 

Table VI using the two-dimensional table lookup routine DISCOT, given T-j-^ 
and the log to the base ten of the ratio of total pressure at Station 3 at 
the compressor face and a reference turbojet discharge region static pressure 
Pj^ which may be the underwing pressure or the freestream static pressure. 

The above pressure ratio is obtained by first determining 
Table IV, given M^, and multiplying by ^ 

The fuel flow rate, maximum thrust and specific fuel consumption are then 
computed as follows: 

^MAX = \ 


Where is the maximum net thrust of the turbojet given by the 

expression 




- Ram Force 

MAX 


Note that the specific fuel consumption given above has units of (Lb . Fuel/Sec) / 
Lb. Thrust. The non-standard time unit is employed to be compatible with 
subsequent weight computations in NSEG. The fuel/air ratio is obtained 

from Table VII, given T-j-^- 

The program is constructed to allow thrust required in the flight direction 
to be determined either by a turbojet throttle setting input or a thrust 
required input. If the throttle setting option is exercised, thrust required 
is computed as 


where Tj is the turbojet throttle setting input. If F|^ is input greater 
than the net thrust which the turbojet can provide in the combined engine 
mode, the excess thrust will be provided by the ramjet, if poss.ible. Prior 
to entering the turbojet iteration loop the maximum turbojet net thrust is 
computed 


^^N^MAX 


F,,^^ - RAM DRAG 
MAX 


and miscellaneous iteration parameters are initialized. 

The purpose of the turbojet iteration loop is to provide the net thrust 
required, if possible, totally by throttling the turbojet. This procedure 
requires that spillage drag, a function of turbojet throttle setting, be 
accounted for in the net thrust, i.e.. 



(Ftot)tj * 


RAM DRAG - SPILLAGE DRAG 


3S 


In the above expression y is the angle between the wing under surface 
and the freestream direction, RAM DRAG is the component of RAM FORCE 
in the freestream direction and SPILLAGE DRAG includes drag components 
due to inlet spillage and excess captured air spillage, if any. 

The logic within this loop first tests the turbojet throttle setting. 

If this value is less than one, iteration commences. Otherwise, turbo- 
jet throttle setting and total thrust are adjusted to their maximum 
values and a check for the combined engine option is made, if the combined 
engine option is not requested, an optional message will indicate that 
the thrust requested is higher than that which can be provided. Within 
the iteration loop, throttled turbojet performance parameters are computed 
by first determining the throttled engine SFC ratio from Table VIII, 
given the throttled net thrust to maximiam net thrust ratio. The throttled 
SFC, fuel flow rate, and airflow rate then become simply 

SFC = (SFC)^ * SFC/(SFC)^ 

= SFC * 

Since the throttled net thrust to maximum net thrust ratio required to 
compute turbojet airflow rate is a function of turbojet ram force which 
itself is a function of turbojet airflow rate, an iteration loop 
encompassing the above computations is required to ensure that the 
turbojet airflow rate used to compute the net thrust ratio is equal to 
that used to compute the ram drag. 

If as computed above is less than the previously computed inlet captured 
airflow rate, the excess is spilled by the turbojet and results in either 
a spillage drag, a ramjet thrust, or both. If W^;^ is greater than the inlet 
captured airflow rate, the excess is assumed to be available from suck-in 
doors . 

If, (1) the turbojet throttle setting equals one and the combined engine 
option is selected or (2j a desired turbojet throttle setting has been 
specified by the subroutine input data the net turbojet thrust available 
in the freestream direction is computed by the expression: 

FTj = CFTOT!)-j. * COSY “ ram DRAG - INLET SPILLAGE DRAG 
J 

and logic flow exits to ramjet calculations. 

If the turbojet throttle setting is less than one, the spilled air lift and 
drag components due to spillage of excess captured air (engine spillage) 
are computed as shown in Figure 8. The available turbojet thrust 
is updated to account for engine spillage drag by the equation 
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= (F„_) * cosy - RAM DRAG - INLET SPILLAGE DRAG - 

TJ TOl xj engine SPILLAGE DRAG 

A maximum throttle setting check followed by a satisfactory thrust level 
check are performed; satisfaction of either results in a return to the 
calling routine. The iteration loop following no exit is a simple 
application of Newton's method which attempts to satisfy the required 
thrust level within one tenth of one per cent. A maximum of ten iterations 
is allowed to complete this requirement. The independent variable in 
this process is the total turbojet thrust. In the event convergence is 
not achieved in ten iterations, an error message is printed and a return 
to the calling program is executed with the engine parameters computed 
on the last iteration. 


6. 2. 1.6 RAMJET AIRFLOW, FUEL FLOW, AND THRUST COMPUTATIONS 

The ramjet computation sequence is constructed similar to that for the 
turbojet. The maximum ramjet air flow rate to the full capture inlet 
air flow rate is computed from Table X given the total temperature at 
Station 2, ^t 2 ' maximum ramjet air flow is then, simply, 


^^a)rj max (Wa)full * ^'^^^full 

The ramjet actual fuel/air ratio divided by the stoichiometric 
fuel/air ratio and ramjet specified impulse are determined from 
Table IX as a function of freestream Mach number. Engine spillage drag, 
if any, ramjet and total fuel flow rates, air flow rates, and net 
available thrust levels are computed in the ramjet iteration loop via 
the equations 


* ^^A^RJ * ^^f/^A^Stoichiometr1c 

stoichiometric 

Fjy = (Wf)pj * Up * cos Y - SPILLAGE DRAG 

(W£)tot = ^ 

(Wa)tot ' ('*A^RJ * 
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Where the ramjet air flow rate above is either due to the air flow 
spilled from the turbojet or the maximum ramjet air flow rate, which 
ever is less. In the latter case, the excess air is spilled at pressure 
?2 and engine spillage forces are computed. The term SPILLAGE DRAG 
above consists of only engine spillage drag in the combined engine option 
but includes inlet spillage drag in the ramjet only option. The ramjet 
iteration loop again applies Newton's method to satisfy, if possible, 
the thrust level requested within .1 per cent using the ramjet air 
flow rate as the independent variable. If this condition cannot be 
satisfied, an appropriate message is printed and a return to the calling 
routine is executed with the engine parameters last computed. A flowchart 
of the combined engine subroutine (ENGINE) appears in Figure 10. 
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6.2.2 Usage 


The turbojet, ramjet and/or combined engine propulsion option (ICASEN = 5) 
provided by the engine simulation in subroutine ENGINE requires the input 
data presented in Section 6.2.3. This data describes "base" turbojet and , 
ramjet engines and relevant airframe geometry. These base engines may be 
scaled, if desired, by utilizing NSEG input option (21) described in Section 
9.21. Scaling must be accomplished in the first mission segment. 

Segment (2). (Segment (1) may never be used as a mission segment; 
the, data array locations which would be employed are utilized for other 
mission parameters.) 

Turbojet scaling of the base engine is accomplished by (1) specifying 
a desired net turbojet thrust (LB) at zero velocity at sea level or 
(2) specifying the sea level corrected airflow rate (LB/sec) for the 
desired engine. The ramjet is scaled by varying the inlet area to 
provide a specified ramjet net thrust at cruise Mach number and altitude. 

A ramjet thrust margin may also be specified- If option (21) is not 
exercised in segment (2) the engines utilized will be described by 
the base engine data and the input inlet area, AINLET. 

Maximum and throttled (specified) power settings can be employed with the 
turbojet, ramjet and/or combined engine propulsion option. If throttling 
is not specified the thrust provided will be such that the thrust component 
in the freestream direction is equal to that requested (equal to vehicle 
drag in the unaccelerated mode). The ENGINE subroutine is called by the 
general propulsion subroutine, ENDAT. An iteration is provided in 
ENDAT which automatically alters the vehicle angle of attack until 
the sum of the vehicle aerodynamic lift and the thrust component in 
the lift direction is equal to the lift necessary for the flight 
condition specified (equal to the vehicle weight less centrifugal lift 
for level flight) . 
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6.2.3 Input Data for Turbojet, Ramjet, and/or 
Combined Engine Option 


In addition to the segment data described in Section 9.21 the following 
data is required to utilize the turbojet, ramjet and/or combined engine 
propulsion option. 


SYMBOL 

DESCRIPTION 

INDE03 

Indicator for usage of engine data table 3. 

INDE04 

Indicator for usage of engine data table 4. 

INDE05 

Indicator for usage of engine data table 5. 

INDE06 

Indicator for usage of engine data table 6. 

INDE07 

Indicator for usage of engine data table 7. 

INDE08 

Indicator for usage of engine data table 8. 

INDE09 

Indicator for usage of engine data table 9. 

INDEIO 

Indicator for usage of engine data table 10. 


NOTE: All the above indicators must be set 


positive for this option. 

INDCLZ 

Indicator for usage of C^q Table. 

INDCLA 

Indicator for usage of Cl Table. 


SYMBOL 

DESCRIPTION 

AINLST 

Engine Inlet Area (initial guess if ramjet is 
to be scaled) (ft^) . 

DELWNG 

Angle from vehicle reference line to wing lower 
surface (degrees) . 

DELIN 

Inlet wedge angle (degrees) . 

AMAXTJ 

Maximum Mach number for turbojet operation. 

AMINRJ 

Minimum Mach number for ramjet operation. 
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SYMBOL 


DESCRIPTION 


IE03X 

IE04X 

IE05X 

IE06X 

IE06Z 

IE07X 

IE08X 

IE09X 

lElOX 

ICLAX 

ICLAY 

ICLZX 

ICLZY 


Number of Mach number (Mq) values in ETAB03. 

Number of Mach number (M^) values in ETAB04. 

Number of Temperature values (Tt^) ETAB05, 

Number of Temperature values (Tt^) ETAB06. 

Number of Pressure ratio values Logjo 
in ETAB06. 

Number of Temperature values ETAB07. 

Number of Throttle setting values in ETAB08. 
Number of Mach number (Mq) values in ETAB09. 
Number of Temperature values ETABIO. 

Number of Mach number values in CLATAB. 

Number of Altitude values in CLATAB. 

Number of Mach number values in CLZTAB. 

Number of Altitude values in CLZTAB. 


IE0P 


Engine 
0 - 


1 - 
2 - 
3 - 


option Selector 

use Turbojet only if M < AMINRJ 

use Combination if AMINRJ < M < AMAXTJ 

use Ramjet only if M >. AMAXTJ 

use Turbojet only 

use Combination 

use Ramjet only 




SYMBOL 


DESCRIPTION 


ETAB03 


ETAB04 


ETAB05 


ETAB06 


Engine Capture Area Schedule 
Ratio = fCMj) 

ETAB03 = mJ,-— R atio^, — , Ratio^ 

N = IE03X = Number of Data pairs 

Inlet Pressure Recovery Schedule 
Ratio = f(Mj) 

ETAB04 = Ratio^ , ---Ratio^ 

N = IE04X 

Base Turbojet Corrected Airflow Schedule (LB/SEC) 
Airflow = f(T-{-^) 

ETAB05 = T^^,---,Tt^, Airflow^,---, Airflow^ 

N = IE05X 

Turbojet Gross Thrust /Airflow rate Schedule 
(LB Thrust/ (LB Airflow/SEC) ) 

FGWA = Logjg (Ptj/PN)) 

ETAB06 = 

fgwaL— ,FGWA^, 

FGWA^, — .FGWaJJ. 

N = IE06X 


M = IE06Z 




SYMBOL 


DESCRIPTION 


ETAB07 


ETAB08 


ETAB09 


ETABIO 


Turbojet Fuel/Air Ratio 
WFWA = fCTtj) 

ETAB07 = Tt^,— ,Tt^, WFWA^ , — -WFWA^ 
3 3 

N = IE07X 



Ramjet Maximum Airflow Rate Schedule 
WOWRAM = f(Ttj) 

ETABIO = Tt , W0WRAM^,-—, WOWRAM^ 

‘"7 -Z 


N = lElOX 



SYMBOL 


DESCRIPTION 


CLZTAB 


CIATAB 


C, at zero angle of attack 

Li 

C = f(Mo,h) 

CLZTAB = mJ, — -,Mg, h^,— -h , 

c 1 — c ** 

C ^ — .Cl“ 

^Om Ow 


N = ICLZX 


M = ICLZY 


Lift Curve Slope Table 
Cj - f(M ,h) 

CLATAB = mJ,— -,M o, h^,- 
C 1 — C ^ 

Xi’ X^’ 

I c ^ 


N = ICLAX 


M = ICLAY 





7. 


AERODYNAMICS 


7 . 1 Introducti on 

Aerodynamics options in NSEG compute the total vehicle drag force given 
the required lift force. Lift force required is determined in the various 
mission segment options available in NSEG , There are two basic aZean 
vehicle aerodynamics options; in each option a total of three aerodynamic 
systems may be defined. Aerodynamic specifications in the two clean vehicle 
options are 

Cd^ = 

and 

CDj = CDojCh, 

j = 1, 2, 3 


* ^Dinducedj «) J j = 1, 2, 3 
M) + Klj ♦ K2jCCL - CL^„(M)f + Cl^; 


In either option the user may optionally sum the system drag components so 
that 

Cd = Cdi + C02 + 

or select the lowest drag configuration 

Cj) = Min [CDj , ^D 2 * ^Dj3 

or finally, any one of the clean vehicle aerodynamic systems may be specified 


= ^Dj i j = 1 or 2 or 3 

Store and tank drags together with their associated pylons may also be 
specified. Three systems of stores are available; each system may involve 
any number of store pairs and store pylon pairs. The aerodynamic form 
employed for store drag is 

^D5 “ I (^DSj ^Sj ^Dsp • ^SPj ) * j = 1 » 2 , 3 

j ^ 

where 


NSj 

^DsPj 

NsPj 


drag of single store pair, type j 

number of store pairs , type j 

drag of single pair of store pylons , type j 

niomber of store pylon pairs , type j , whose drag is not included 
in the clean vehicle deag tables 
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Similarly, tank system drags are written in the form 



where the notation is similar to that employed for stores. 

7.2 Aerodynamic Subroutines 

The subroutines of particular significance in aerodynamic drag calculations 
are 

DRAG - The control program for aerodynamic drag calculations 

SACNSl - First clean vehicle aerodynamics option available in 
NSEG 

SACNS2 - Second clean vehicle aerodynamics option available 
in NSEG 

SACSP - Store and store pylon drag calculation 
SACTP - Tank and tank pylon drag calculation 
Flow charts for these five subroutines are presented in Figures 10 through 14. 


7 . 3 Data Input 
7.3.1 General 


SYMBOL DESCRIPTION 

INDAER Aerodynamic option indicator: 

= 8, Option 8, tabular 

= 9 , Option 9 , derivatives 

^ 8 or 9 , error 

CASE Aerodynamic c^lse indicator: 

= -1, Use smallest drag from Systems 1, 2, and 3 
= 1, Use System 1 

= 2, Use System 2 

= 3, Use System 3 

= 4, Sum System 1, 2, and 3 
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ATAB40 

ATAB41 

ATAB42 

ATAB43 

ATAB44 

ATAB45 

IA40X 

IA41X 

IA42X 

IA43X 

IA44X 

1A45X 

IA40Y 

IA41Y 

IA42Y 

IA43Y 

IA44Y 

IA45Y 

INDA40 

INDA41 

INDA42 

INDA43 

1NDA44 

INDA45 


7.3.2 Aerodynamics Option 8 


Zero lift drag, System 1; = Cog^Ch, M) 

Zero lift drag. System 2; CDq 2 = K) 

Zero lift drag. System 3; CDqj = COq^CH, M) 

Drag due to lift. System 1; M) 

Drag due to lift. System 2; = CDj^ 2 (CL» 

Drag due to lift. System 3; M) 


Number 

Number 

Number 

Number 

Number 

Number 

Number 

Number 

Number 

Number 

Number 

Number 


of altitude value in ATAB40 
of altitude value in ATAB41 
of altitude value in ATAB42 
of value in ATAB43 
of Cl value in ATAB44 
of Cl value in ATAB45 
of Mach number values in ATAB40 
of Mach number values in ATAB41 
of Mach number values in ATAB42 
of Mach number values in ATAB43 
of Mach number values in ATAB44 
of Mach number values in ATAB45 


Table control indicator ATAB40 
Table control indicator ATAB41 
Table control indicator ATAB42 
Table control indicator ATAB43 
Table control indicator ATAB44 
Table control indicator ATAB45 



7,3.3 Aerodynamics Option 9 


ATAB46 

ATAB47 

ATAB48 

ATAB49 

ATAB50 

ATAB51 

ATAB52 

ATAB53 

ATAB54 

ATABSS 

ATAB56 

ATAB57 

ATABSB 

ATAB59 

ATAB60 

IA46X 

IA47X 

IA48X 

IA46Y 

IA47Y 

IA48Y 

INDA46 

INDA47 


Zero lift drag. System 1, Cdqj^ = M) 

Zero lift drag. System 2, Cdq 2 = CDQ2Ch> M) 

Zero lift drag. System 3, M) 

Coefficient of , System 1. I'll - Kii(M 3 
Coefficient of Cj^ , System 2, K12 = Ki2(M) 

Coefficient of Cl , System 3, K13 = Ki3(M) 

Coefficient of (Cl - * System 1, K21 == K2 l(M) 

Coefficient of (Cl - > System 2, K22 = ^22(^3 

Coefficient of (Cl - Cj^^) , System 3, K23 = K23(M) 

Coefficient of C^ , System 1, K32 = K3i(M) 

Coefficient of Cl » System 2, K32 - K32(M) 

Coefficient of Cl » System 3, K33 = K33(M) 

Minimum drag Cl, System 1, = ^Lniin.l^^^ 

Minimum drag Cl, System 2, 2 ^Ljnin.2^^^ 

Minimum drag Cl, System 3, 3 = 3(M) 

Number of altitude values in ATAB46 
Number of altitude values in ATAB47 
Number of altitude values in ATAB48 
Number of Mach number values in ATAB46 
Number of Mach number values in ATAB47 
Number of Mach number values in ATAB48 

Table control indicator ATAB46 
Table control indicator ATAB47 
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INDA48 

Table 

control 

indicator 

ATAB48 

INDA49 

Table 

control 

indicator 

ATAB49 

INDASO 

Table 

control 

indicator 

ATAB50 

INDA51 

Table 

control 

indicator 

ATAB51 

INDA52 

Table: 

control 

indicator 

ATAB52 

INDA53 

Table 

control 

indicator 

ATAB53 

INDA54 

Table 

control 

indicator 

ATAB54 

INDA55 

Table 

control 

indicator 

ATAB55 

INDA56 

Table 

control 

indicator 

ATAB56 

INDA57 

Table 

control 

indicator 

ATABS7 

INDA58 

Table 

control 

indicator 

ATAB58 

INDA59 

Table 

control 

indicator 

ATAB59 

INDA60 

Table 

control 

indicator 

ATAB60 



8. LANDING AND TAKE-OFF CALCULATION 


NSEG contains an approximate landing and take-off calculation. This model 

has its own aerodynamic (DATCOM) and propulsion data. It is not a mission 

segment option at the present time; however, it is used to perform landing 

and take-off analyses independently of the mission sequence. The landing 

and take-off analysis options are selected by data input via Option 17. 

Where detailed high lift aerodynamic values are known, the DATCOM 

calculations may be bypassed. In this case the values of Cl , ^Lmax* 

and Ow.v (CLAT0T, CLMXT0, ALMXT0) are specified directly in ^ata input. 

MAX 


wo 


AG 

AMAXG 

ARTHEO 

BFOBEX 

CDO 

CDPLG 

CFOCAV 

CLO 

CLATOT 

DFK 

FLAPDT 

FNO 

FFR 

GRNDFT 

SFC 

SPLG 


8.1 Data Input, 

Take-Off 

Weight at start of take-off, in pounds. NOTE that if 

WO < 0.0 

take-off calculation is omitted. 

Ground roll body angle, degrees 

Maximum ground rotation angle, degrees 

Theoretical wing aspect ratio 

Flap span to exposed wing span ratio 

Drag coefficient at zero lift 

Flat plate drag coefficient of gear 

Average flap chord to wing chord ratio 

Lift coefficient at zero angle of attack 

Total linear lift curve slope, per degree 

Induced drag factor 

Take-off flap deflection, per degree 

Take-off thrust, in pounds 

Take-off fuel flow rate, in pounds per hour 

Take-off ground rolJ friction coefficient 

Internal, specific fuel consumption, (hour 

Flat plate area of gear 
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SREF 


TRATIO 

WO 

WINCID 

WSLEFC 

XISO 

The landing 
WL 

CDPCHT 

FLAPPDL 

FNL 

GRNDFL 

SPCHUT 

WL 

XMAXG 


Wing reference area 
Wing taper ratio 
See above 

Wing incidence, degrees 

Wing leading edge sweep angle, degrees 

Internal, specific impulse at take-off, in seconds 

8.2 Data Input, Landing 

calculation uses the take-off data with the following modifications: 

Landing weight. NOTE that if WL ^ 0.0, no landing perfor- 
mance calculation is performed. 

Flat plate drag coefficient of landing chute 

Landing flap deflection 

Landing thrust 

Landing ground roll friction coefficient 
Flat plate area of landing chute 
See above 

Maximum ground roll distance, in feet. It is ignored if 
equal to 0.0; otherwise the thrust required to achieve the 
desired ground roll is computed. 
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9. MISSION SEGMENT OPTIONS 


NSEG contains twenty-one mission segment options listed below. 

1. Climb or descent, includes approximately optimal paths 

2. Maximum rate of climb or maximum range descent 

3. Radius mission 

4. Cruise-climb to specified weight 

5. Cruise-climb for specified range or time 

6. Constant altitude cruise to a specified end weight 

7. Constant altitude cruise for a specified range 

8. Constant altitude cruise for a specified time 

9. Buddy refuel 

10. Mach-altitude-weight transfer 

11. NSEG logic option 

12. Instantaneous weight change 

13. Instantaneous Mach-altitude change 

14. General purpose and point condition calculations 

15. Iterate to fly a specified total distance 

16. Climb or accelerate 

17. Fuel weight change 

18. Not used 

19. Not used 

20. Fuel allowance 

21. Turbojet, ramjet and/or combined engine scaling option 

22. Range matching generalized iterative control option 

Data input descriptions and typical cases are described in the remainder of 
this section. 
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9.1. Option 1. Climb or Descent Path 


This option calculates a climb or descent path between two Mach-altitude points, 

(Hq, AMS) -► (HE, AME) 

where 

(Hq, AMS) - the initial Mach-altitude point and is taken as the current Mach 
number and altitude point 

(HE, AME) - final Mach-altitude point 

If the value of the final Mach-altitude point is known, it can be specified 
directly as data. In some cases the final Mach-altitude point is the result 
of a calculation in a previous segment. In this case. Option 10 would be 
used to transfer the Mach number and altitude from the segment in which it was 
calculated into the segment in which the climb-descent path is to be computed. 
Range, time, and fuel used in the climb or descent are accounted for when 
Option 1 is specified. It also allows either forward or reverse flight path 
calculation. 

When using Option 1, the intermediate Mach-altitude points may be selected in 
any one of seventeen ways. Figure 16 shows the various methods available for 
selecting the intermediate Mach-altitude points along the flight path. If this 
option is used for climb paths, thrust equal to maximum should be specified; 
for descent paths, thrust equal to minimum should be specified. The flight 
path is adjusted at the beginning and end to meet the specified end points. 


Data for Option 1 . 

The following is a list of all possible data statements for Option 1. The 
information in ■itdli-cs is optional . Note that v and n are decimal numbers and 
must contain a decimal point. 

end MACH number =v, 

end ALTITUDE=n, (feet) 

USE MAX EDOT ae baeis for aeteotton of Maoh-^ attitude points 

USE MAX EDOT/MDOT as basis for seleoticn of Maoh- attitude points 

USE MAX L/D aa basis for setection of Maoh~attitude points 

USE MAX RANGE FACTOR as basis for setection of Mach~attitude points 

USE MAX L/(FN-DRT) as basis for setection of Mach^attitude points 
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sectrdh for Machr- altitude point using CONSTANT ENERGY 
eearoh for Maah-altitude point using CONSTANT ALTITUDE 
search for Maah-altitude point using CONSTANT MACH NUMBER 
fly CONSTANT Q path 
fly LINEAR Maah-altitude path 

REVERSE calculation (with respect to flight path) 

* 

forward calculation (with respect to flight path) 


Example Data Set for Option 1 . 


END SEGMENT (3) 


END MACH NUMBER = 2.00, 


END ALTITUDE * 45000.0, 


USE MAX EDOT 


CONSTANT ENERGY 


THRUST = MAXIMUM 


SEGMENT(3), OPTION(l) 

The maximum E path flown by this option is shown in Figure 16. 

A list of all flight paths available in Option 1 is presented in 
the following figure. 



* An entire statement appearing in italics denotes the default option. 
The statement has no effect on the program. 
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9.2. Option 2. Maximum Rate of Climb or Maximum Range Descent Path 

Option 2 calculates a^ climb or descent such that the vehicle follows a flight 
path for maximum rate of climb or minimum rate of sink. The program determines 
the speed for maximum rate of climb at intermediate altitudes. The initial 
condition is the Mach- altitude- weight at the end of the previous segment. 

The flight path is adjusted at the beginning and at the end to meet the 
specified end points. 


Data for Option 2 . 

The following is a list of data statements for Option 2. The information in 
italics is optional . Note that v and n are decimal numbers and must contain 
a decimal point. 

end MACH nwrbev =v, 

end ALTITUDE=*n, (feet) 

catautate DESCENT path 

catoulate cUrrb path * 

REVERSE calculation (with respect to flight path) 
foicward calculation (with respect to flight path) * 

Example Data Set for Option 2 . (Climb Path Calculation) 


/ 



END SEGMENT (5) 


THRUST=MAX 


END MACH NUMBER 0.76, 


END ALTITUDE = 37500.0, (FEET) 


SEGMENT (5) , OPTION (2) 


♦ An entire statement appearing in italics denotes the default option. 
The statement has no effect on the program. 
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It should be noted that Option 2 was part of the original NSEG program 
of Reference 1 and was included in the current version of NSEG without 
modification in order to retain the same option number as that of the 
original program. 

It is recommended that the maximum E climb of Option 1 be used in place 
of Option 2. 
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9.3. Option 3. Radius Mission Option 


When Option 3 is specified, the program performs an iteration process such 
that the range gained in some sequential portion of a mission is made equal 
to the range gained in some other sequential portion of the mission. The 
range adjustable segment must be either an Option 5 segment or an Option 7 
segment. 


Data for Option 3 . 

The following is a list of data statements for Option 3; information in 
italics is optional. Note that i, j, k, t, and m must be integer numbers; 
therefore, decimal points should not be used. 

START fivat LEG at aegmentii) 

END firat LEG at aegment (j) 

START RETURN LEG at aegment (k) 

END RETURN LEG at aegment (A,) 

ADJUST range in aegment (m) 


Example Data Set for Option 1 . 



END SEGMENT (30) 


ADJUST RANGE IN SEGMENT (8) 


END RETURN LEG AT SEGMENT (29) 


START RETURN LEG AT SEGMENT (17) 


END FIRST LEG AT SEGMENT (12) 


START FIRST LEG AT SEGMENT (3) 


SEGMENT (30), OPTION (3) 
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9.4. Option 4. Cruise Climb to a Specified Weight 


Mission option 4 calls for the computation of a cruise climb flight segment 
from one weight to another at aonstccnt Mach nwrher and constant It ft 
ooefftaient. This type of cruise climb is also called a constant w/6 
cruise since the vehicle weight divided by the pressure ratio is a constant. 
The flight path for this option may be computed using one step or five 
steps . 

Data for Option 4. 

The following is a list of data statements for Option 4; information in 
italics is optional. Note that v is a decimal number and therefore must 
contain a decimal point. 

end WEIGHT=v, pounds 

make five FLIGHT PATH corrections 

DO NOT make FLIGHT PATH corrections 


Example Data Set for Option 4 . 


/ 

/■ 

/ 

/■ 

/ 


END SEGMENT (15) 


THRUST = THRUST REQUIRED 


MAKE FIVE FLIGHT PATH CORRECTIONS 


END WEIGtfT = 335000.0, POUNDS 


SEGMENT (15), OPTION (4) 
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9.5. Option 5. Cruise Climb for Specified Range or Time 


Mission Option 5 calls for the computation of a cruise climb flight segment 
for a specified range or a specified time. As in Option 4 this cruise 
climb flight path is at constant Mach and constant lift coefficient. 


Data for Option 5. •" 

The following is a list of data statements for Option 5; information in ■ 
italics is optional. Note that v is a decimal number and must contain a’ 
decimal point. V 


cruise TIME-v, (hours) 

cruise RANGE=v, (n» mites in this segment) 

TOTAL RANGE fiown including this segment =v, (n,m*) 
NO RANGE credit 

REVERSE calculation (with respect to flight path) 
forward calculation with respect to flight path * 


Example Data Set for Option 5 . 


/ 


END SEGMENT (17) 


THRUST = THRUST REQUIRED 


TOTAL RANGE FLOWN INCLUDING THIS SEGMENT = 1500.0, 


SEGMENT (17), OPTION (5) 


(N.M.) 


* An entire statement appearing in italics denotes the default option. 
The entire statement has no effect on the program. 
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9.6. Option 6. Constant Altitude Cruise to a Specified End Weight 


Option 6 calculates a constant altitude cruise path to a specified end weight. 
The constant altitude cruise may be performed at a constant Mach number or 
at a constant lift coefficient. The initial condition is the Mach-altitude- 
weight at the end of the previous segment. 

This option provides the capability of flying the path backwards by 
specifying a weight that is larger than the vehicle weight at the end 
of the previous segment. 


Data for Option 6 . 

The following is a list of data statements for Option 6; information in italics 
is optional. Note that ri is a decimal number and must contain a decimal 
point. 

WEIGHT at end of oruiae s*n> POUNDS 
oruise at constant LIFT ooeffiaient 
cruise at constant Mach number 


Example Data Set for Option 6 , 


END SEGMENT (21) 


THRUST = THRUST REQUIRED 


WEIGHT AT END OF CRUISE = 76000.0, POUNDS 


CRUISE AT CONSTANT MACH NUMBER 


SEGMENT (21), OPTION (6) 


* An entire statement appearing in italics denotes the default option. 
The statement has no effect on the program. 
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9.7. Option 7. Constant Altitude Cruise for a Specified Range 

This mission option calls for the confutation of a constant altitude cruise 
flight segment for a specified range at constant Mach number. 

Data for Option 7 . 

The following is a list of data statements for Option 7; information in 
italics is optional. Note that v is a decimal number and must contain a 
decimal point. 


RANGE flcwn in this segment -v, (n,m.) 

TOTAL RANGE flown including this segment =v, (n.m, ) 
REVERSE calculation (with ipespect to flight path) 
foxward calculation (with respect to flight path)* 


Example Data Set for Option 7 . 



END SEGMENT (29) 


TOTAL RANGE = 2700.0, (N.M.) 


THRUST = THRUST REQUIRED 


SEGMENT (29), OPTION (7) 


* An entire statement appearing in italics denotes the default option. 
The statement has no effect on the program. 
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9.8. Option 8. Constant Altitude Cruise for a Specified Time 


This mission option calls for the conq)utation of a constant altitude cruise 
flight segment for a specified time at either constant Mach number or 
constant lift coefficient. This option may also be used to compute fuel 
required for loiter by specification of no range credit in the data set. 


Data for Option 8 . 

The following is a list of data statements for Option 8; information in 
italics is optional. Note that v is a decimal number and must contain a 
decimal point. 


oruise TIME=v, (hours) 

NO RANGE oredit 

cruise at constant LIFT coefficient 
cruise at constant Mach number * 

REVERSE calculation (with respect to flight path) 
forward calculation (with respect to flight path) 


Example Data Set for Option 8 . 



END SEGMENT (24) 


CRUISE TIME =2.5, (HOURS) 


CONSTANT LIFT COEFFICIENT 

i 


THRUST = THRUST REQUIRED 


SEGMENT (24), OPTION (8) 


* An entire statement appearing in italics denotes the default option. 
The statement has no effect on the program. 
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9.9. Option 9. Buddy Refuel 


Option 9 will determine where the ''optimum'* in flight buddy refuel should 
take place and how much fuel will be transferred. The user must input 
three points of the tanker off load curve in the DATAl data block. The 
names of the tanker off load data are shown below. 


POINT 

(RANGE (N.M.) 

FUEL AVAILABLE (POUNDS) 

1 

RANGE 1 

FUEL 1 

2 

RANGE 2 

FUEL 2 

3 

RANGE 3 

FUEL 3 


The tanker off load data will be curve fitted to determine the coefficients 
of a quadratic equation and the quadratic equation will then be used in the 
mission option to determine fuel available as a function of range. 

The cruise out to the refuel point will be one of the following: 


1. 

Constant 

Mach, 

constant 

Cl 

2. 

Constant 

Mach, 

constant 

altitude 

3. 

Constant 

Cj^, constant altitude 


Option 9 is similar to Option 4 and Option 6 except that the segment end 
weight will be determined by the refuel criteria. The vehicle will cruise 
out iintil one of the following conditions is satisfied: 

1. The fuel receivable equals the fuel available. 

2. The distance flown equals the maximum allowable 
distance to refuel if such a constraint exists 

3. The minimum in-flight weight of the receiver is 
reached. 

If the fuel receivable exceeds the fuel available at the start of the mission 
segment in which Option 9 is used, transfer of the fuel will take place at 
that point in the mission. 


Data for Option 9 . 

The following is a list of data statements for Option 9; information in 
italics is optional. Note that v and q are decimal numbers and must contain 
a decimal point. 
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maximum RANGE to refuel =t>, (n.m,) 
no maximum refuel range is to be considered * 
cruise at CONSTANT CL, constant altitude 
cruise at CONSTANT MACH, constant CL 
CRUISE AT CONSTANT MACH, constant altitude 
imximum in flight WEIGHT=v, (pounds) 
fuel transfer RATE=n, (pounds /minute) 

Note that if the fuel transfer rate is not specified or is specified as a 

value less than or equal to zero, the fuel transfer will be instantaneous. 


Example Data Set for Option 9 . 



END SEGMENT (21) 


MAXIMUM RANGE TO REFUEL =* 2200.0, (N.M.) 


CRUISE AT CONSTANT MACH, CONSTANT CL 


MAXIMUM IN FLIGHT WEIGHT = 450000.0, POUNDS 


SEGMENT (21), OPTION (9) 


* 


An entire statement appearing in italics denotes the default option. The 
statement has no effect on the program. 



9.10. Option 10. Mach-Altitude-Weight Transfer 


Option 10 transfers values calculated in any segment into some other segment; 
for example, the end Mach-altitude point for a climb may be the result of an 
Option 14 point calculation. Option 10 would be used to transfer the Mach- 
altitude values from the segment in which the information was calculated into 
the segment in which the climb is to be performed. 


Data for Option 10 . 

The following data statement is used for Option 10. Note that v and n must 
be integer numbers and cannot contain a decimal point. 


/MACH 

TRANSFER {ALTITUDE 

/weight 


FROM SEGMENT (v) TO SEGMENT (n) 


Example Data Sets for Option 10 . 



END SEGMENT(17) 

TRANSFER MACH ALTITUDE WEIGHT FROM SEGMENT (3) TO SEGMENT (30) 
SEGMENT(17), OPTION(IO) 

• • • NOTE MACH ALTITUDE AND WEIGHT MAY BE IN ANY ORDER 


END SEGMENTC16) 

TRANSFER ALTITUDE FROM SEGMENT (9) TO SEGMENT (21) 
SEGM£NT(16), OPTION(IO) 

• .. NOTE ONLY VALUES SPECIFIED WILL BE TRANSFERRED 


END SEGMENT(15) 

TRANSFER ALTITUDE MACH FROM SEGMENT(6) TO SEGMENT(19) 
SEGMENT(IS), OPTION(IO) 

THREE EXAMPLES OF USING OPTION(IO) 
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9.11 Option 11. NSEG Logic Option 


This option allows for choosing one of two sets of flight conditions 
(Mach, altitude, and weight) by deciding which has the largest or 
smallest value of the variable used as the basis of selection. The 
user specifies the Mach, altitude, and weight for flight condition 2 
and the segment in which flight condition 1 was calculated. The user 
also specifies the basis of selection: 

Maximum Mach 
Minimum Mach 
Maximum Altitude 
Minimum Altitude 
Maximum Weight 
Minimum Weight 


Data for Option 11 . 

The following is a list of data statements for Option 11; information in 
italics is optional. Note that v, n, and X are decimal numbers and must 
contain a decimal point; A is an integer number and cannot contain a 
decimal point. 

MACH number for flight condition 2 =v, 

ALTITUDE for flight condition 2 (feet) 

WEIGHT for flight condition 2 =X, (-pounds) 

get data for flight condition 1 from SEGMENT (A) 

^ I MAXIMUM/ ( MACH 

CHOOSE on the ba8^8 of j minimum! J ALTITUDE 

' * WEIGHT 

Example Data Set for Option 11 . 


END SEGMENT 


( 6 ) 



GET FLIGHT CONDITION 1 FROM SEGMENT 


CHOOSE ON MAXIMUM MACH 


(4) 


MACH = 2.25, 


/ WEIGHT = 78000.0, (POUNDS) 


/ 


ALTITUDE = 64000.0, (FEET) 
SEGMENT (6), OPTION (11) 
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9.12. Option 12. Instant . Weight Change 


This option allows instant changes in the vehicle weight; for example, 
changes that would occur with the dropping of an internal payload. 


Data for Option 12 . 

The following data statement is used for Option 12. Note that v must be 
decimal point. The weight of the vehicle will be computed by subtracting 
V from the current vehicle weight. 


WEIGHT change =v, (POUNDS) 


Example Data Set for Option 12 . 


A 

A 

A 

A 

A 

A 

A 


END SEGMENT (9) 


WEIGHT CHANGE = 10000.0, POUNDS 


••• DROP 10000.0 POUND WEAPON 


SEGMENT (9), OPTION 12 


END SEGMENT (8) 


WEIGHT CHANGE = -74000.0, POUNDS 


. .. INSTANT REFUEL OF 74000.0 POUNDS 


SEGMENT (8), OPTION (12) 



9.13. Option 13. Instant Mach<>Altitude Change 


This option provides the user with a method of changing the Mach-altitude- 
weight point of the vehicle. Option 13 changes any combination of Mach 
number, weight, and altitude. 


Data for Option 15 . 

The following is a list of data statements for Option 13; information in 
italics is optional. Note that v and n ^re decimal numbers and must 
contain a decimal point. 

neu MACH number =v, 

new ALTITUDE * Y, (feet) 

new vehicle WEIGHT=n, (pounds) 


Example Data Set for Option 13 . 



END SEGMENT (31) 


NEW MACH NUMBER = 0.42, 


NEW WEIGHT = 38000.0, (POUNDS) 


SEGMENT (31), OPTION (13) 
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9.14. Option 14. General Purpose and Point 
Condition Calculation Option 


Option 14 is used for either of two basic functions: 

1. Restoring of the vehicle to the Mach, altitude, and 
weight condition at the end of the specified segment. 

Vehicle return weight is restored if WEIGHT is set to zero. 

2. Calculating the point condition 

When using the point condition calculation feature of Option 14, the flight 
condition for the calculation may either be specified or taken as the 
vehicle flight conditions at the end of the previous segment. 

The point calculations that can be performed using Option 14 are 

1. Compute best cruise altitude for specified weight eind 
Mach number. 

2. Compute ceiling for the specified rate of climb, Mach 
number, and weight. 

3. Compute Mach number for best lift to drag ratio for 
the specified altitude and weight. 

4. Compute Mach number for lift equal to weight at the 
specified lift coefficient, weight and altitude. 

5. Compute maximum endurance Mach for the specified 
altitude, weight, and maximum lift coefficient. 

6. Compute maximum Mach number for the specified altitude 
and weight . 

7. Compute Mach number for maximum rate of climb at the 
specified altitude and weight. 

8. Compute altitude and Mach number for maximum range factor 
at the specified weight. 

9. Compute Mach number for maximum range factor at the 
specified altitude and weight. 

10. . Confute the various energy maneuverability parameters 
at the specified G level, Mach number, weight, and 
altitude . 
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Data for Option 14 . 


Hollowing is a list o£ data statements for Option 14; information in 
italics is optional. Note that 6, v, n> 3, and X are decimal numbers 
and must contain a decimal point ;£ is an integer number and must not 
contain a decimal point. 

RESTORE vehicle to flight conditions at end of segment ( 1) . 
use flight conditions from previous segment * 

GIVEN MACH number =v, (pounds) 

GIVEN WEIGHT=X, (pounds) . I NOTE: If weight = 0.0, then weight will be 

internally set to return weight. 

GIVEN ALTITUDE^n, (feet) 

GIVEN GLEVEL=6, 

COMPUTE best CRUISE altitude 

COMPUTE CEILING for rate of clirrib =B, (feet/min) 

COMPUTE MACH for best LIFT/DRAG ratio 
COMPUTE MACH for CL=6 , 

COMPUTE MAXirmm endurance MACH for CL=6 , 

COMPUTE rnximum MACH nurrber 
COMPUTE MACH nimher for WJjimum RATE of climb 
COMPUTE MACH, ALTITUDE point for W^^irmm range factor 
COMPUTE MACH for MAXimz^ range factor 

COMPUTE MACH (v) , ALTITUDE point for hlAXimum range factor 

COMPUTE MACH(v), for MAXimum range factor 

COMPUTE ENERGY maneuverability parameters 

Note that for maximum range factor calculation the user can specify the upper 
limit on Mach number; i.e., MACH (v) where v is upper Mach number limit. If 
it is not specified AMLIM is used. 


* An entire statement appearing in italics denotes the default option. The 
statement has no effect on the program. 
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END SEGMENT (18) 


GIVEN GLEVEL =1.5, ALTITUDE = 18000.0, (FEET) 


GIVEN MACH =0.71 


COMPUTE ENERGY MANEUVERABILITY PARAMETERS 


SEGMENT(18), OPTION(14) 




END SEGMENT(17) 


GIVEN ALTITUDE = 18000.0, (FEET) 


GIVEN WEIGHT = 70000.0, (POUNDS) 
COMPUTE MACH FOR BEST LIFT DRAG/RATIO 


SE6MENT(17), 0PTI0N(14) 


/ 

/ 


END SEGMENT(16) 


COMPUTE BEST CRUISE ALTITUDE 


USE FLIGHT CONDITIONS FROM PREVIOUS SEGMENT AND 


SEGMENT(16), 0PTI0N(14) 


/ 


SAMPLE DATA SETS FOR OPTION (14) 



Example Data Set for Option 14 . 



END SEGMENT(21) 


RESTORE VEHICLE TO FLIGHT CONDITIONS AT END OF SEGMENT (3) 


SEGMENT(21), 0PTI0N(14) 


END SEGMENT(20) 


COMPUTE CEILING FOR RATE OF CLIMB =* 100.0, (FEET/MIN) 


GIVEN MACH = 0.72, WEIGHT = 220000.0, (POUNDS) 


SEGMENT(20), OPTION(14) 


END SEGMENT(19) 


COMPUTE MAX ENDURANCE MACH FOR CL 0.7 


GIVEN ALTITUDE = 32000.0, WEIGHT = 270000.0, (POUNDS) 


SEGNIENT(19), 0PTI0N(14) 


( oontinued) 
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9.15. Option 15. Iterate to Fly a Specified Distance 


This option provides an iteration process that forces the vehicle to have 
accumulated a specified distance at the end of this segment. The range 
flown in a specified segment is automatically adjusted by the program until 
the total range flown is within a tolerance of one nautical mile of the 
range specified. 


Data for Option 15. 

The following data statement is used for Option 15; information in italics 
is optional. Note that v is a decimal number and must contain a decimal 
point; £ is an integer number and must not contain a decimal point. 

adjust range in segment (£) to give a total range n,m. 


Example Data Set for Option 15 


END SEGMENT C23) 


ADJUST SEGMENT(IS) SO TOTAL RANGE=1650 .0, N.M. 


SEGMENT(23), OPTION(IS) 



9.16. Option 16. Climb or Accelerate 


With Option 16 two Mach-altitude conditions are specified by the segment 
numbers in which they were calculated. Option 16 makes decisions as per 
MIL-C-Rules to determine if fuel is to be *'bumed” while going from ' 
Condition 1 to Condition 2. 

The decision process is 

If ^2 greater than 

or 

If H 2 equals 

and 

is greater than M, 
then 

The vehicle will fly (burning fuel) to Condition 2 from Condition 1 

If neither of these conditions eixsts, the program will move the vehicle to 
Condition 2 with no change in weight or performance. 

Note that Option 16 cannot be used to descend or decelerate. 


Data for Option 16. 

The following is a list of data statements for Option 16; information in 
italics is optional. Note that I and j are integer numbers and must not 
contain a decimal point. 


obtain condition ONE from segment 
obtain condition TWO from segment 
NO RANGE credit for this segment 
RANGE credit for this segment 
forward catcutation (raith respect 
REVERSE calculation (with respect 


(A) 

(j) 


to flight path) 
to flight path) 


* 


* An entire statement appearing in italics denotes the default <^ti(m. 
The statement has no effect on the program. 
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Example Data. Set for Option 16 . 



END SEGMENT (12) 
RANGE CREDIT 


FORWARD CALCULATION 

GET CONDITION TWO FROM SEQUENT (4) 

GET CONDITION ONE FROM SEGMENT(7) 


SEGMENTC12), 0PTI0N(16) 



9.17. Option 17. Fuel Weight Change for Time, Turn, Combat, etc. 


The fuel weight change may be computed or input depending on the needs of 
the program. If the option of computing the fuel weight change is desired, 
fuel used during warm-up, take-off, loiter, or combat can be computed by 
the program. Landing calculations are also controlled by Option 17. 


Data for Option 17 . 

The following data statements are used by Option 17; information in italics 
is optional. Note that n and X are decimal numbers and must contain a 
decimal point. 

fueZ WEIGHT change (pounds) 

WARM-UP TIME=n, (minutes) 

loiter TIME=n, (minutes) 

n DEG X G TURNS 

COMPUTE COMBAT fuel allowance 

COMPUTE WARM-UP and TAKE-OFF fuel allowance 

COMPUTE loiter fuel allowance 
COMPUTE LANDING performance ' 

Example Data Sets for Option 17 . 


END SEGMENT (27) 


720.0 DEG 4.2 G TURNS 


COMPUTE COMBAT FUEL ALLOWANCE 


THRUST=MAX 


SEGMENT(27), OPTION(17) 


END SEGMENT (2 6) 


LOITER TIME=15.0, (MINUTES) 


COMPUTE LOITER FUEL ALLOWANCE 


THRUST=THRUST REQUIRED 


SEGMENT(26), 0PTI0N(17) 



9.18. Option 18. 


Option 18 is not used in the current version of the NSEG program. 
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I 


9.19. Option 19. 


Option 19 is not used in the current version of the NSEG 


program. 
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9.20. Option 20. Fuel Allowance 

Option 20 computes a fuel allowance for (a) a specified time at a specified 
power setting and (b) a specified tine at a specified thrust to weight ratio. 

Data for Option 20. 

The following data statements are used by Option 20; information in italics 
is optional. Note that n and v are decimal numbers and must contain a 
decimal point. 

GIVEN TIME=v, (minutee) 

confute fuel atlcwanoe for a POWER setting =n, 
compute fuel allowance for tkruat/weight RATI0=n* 


Example Data Set for Option 20. 



END SEGMENT (36) 


POWER SETTING=0.45, 


GIVEN TIME=13.0, (MINUTES) 


SEGMENT(36), OPTION(20) 


END SEGMENT(35) 


THRUST/WEIGHT RATI0=.25, 


GIVEN TIME=17.0, (MINUTES) 


SEGMENT(35), OPTION (20) 
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9.21 Option 21. Turbojet, Ramjet and/or Combined 

Engine Scaling 


Option 21 performs engine scaling computations by performing appropriate 
calls to subroutine ENGINS. Turbojet scaling is performed by computing 
a scaling factor by which the base engine airflow schedule (input) is 
multiplied throughout the flight regime. This factor is computed based 
on a desired net thrust at sea level or a desired turbojet corrected airflow 
at sea level. For an aircraft utilizing turbojet engines only, the cruise 
Mach number and altitude must also be input to allow sizing of the turbojet 
inlet area to meet airflow requirements at the specified cruise condition. 

The ramjet scaling, which sizes the combined engine inlet area, is performed 
at the specified cruise Mach number and altitude to provide net thrust 
equal to drag plus a specified thrust margin at cruise. 

The lift coefficient utilized in drag coefficient calculations is the 
specified cruise lift coefficient (CLC) input in the general data section. 

If CLC is not input (or is input as zero) the lift coefficient for scaling 
is computed using the cruise Mach number and altitude and lift equal to 75% 
of the gross takeoff weight. 

A final option for ramjet scaling is the same as that just mentioned except 
current weight, Mach number, and altitude are used to compute the lift 
coefficient. This option allows the combined engine inlet area to be 
scaled during the mission, if desired. An additional input to this option 
allows the usage of turbojet engine reference pressure (Pj^) to be specified 
as the Wing Undersurface Pressure; if reference pressure is not specified 
ambient static pressure is used. Finally, a combined engine option debug 
print can be requested for engine performance evaluation during climb 
segments. Option 21 must be exercised in Segment (2) if the combined 
engine option is to be employed and engine scaling is desired. If no scaling 
is desired the base turbojet engine airflow schedule and combined engine 
inlet area specified in the general data section will be employed. 


Data for Option 21 

The following data statements are used by Option 21; information in 
italics is optional. Note that all numbers (a, 3, etc.) are decimal 
numbers and must contain a decimal point 


PROPULSIVE mode = DUAL 

SCALE TURBOjet THRUST = a, 
or 


NOTE: If ot <3.0, then a must be the thrust/ 

weight ratio instead of thrust. 


SCALE TURBOjet oovreoted AIRFLOW ^ 3, 


SCALE RA^^'e^ = vamjet thrust MARGIN = e, crutse MACH = X) 
orui-se ALTITUDE = 6 
engine DEBUG print requested 

turbojet reference PRESSURE = underwing pressure 
compute aiphoL for scaling based on current WEIGHT 
END 
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Sample Data Set for Option 21 


SEGMENT (2), OPTION (21) 
'PROPULSIVE MODE = DUAL 


SCALE TURBOJET - SEA LEVEL CORRECTED AIRFLOW = 398.65, (LBS/SEC) 



/CRUISE MACH = 6., ALTITUDE = 120000., (FT) 



9.22 Option 22. Generalized Iteration Control Option 


Option 22 is a generalized iteration loop used to vary a specified parameter 
(independent variable) in a specified mission segment such that a specified 
parameter (dependent variable) either satisfies a given input value to 
within 0.1% or the maximum number of iterations input is exceeded. As 
currently mechanized the dependent variable must be either Mach number, total 
range, or weight. The independent variable must be either Mach number, weight 
or range as specified in the segment where the iteration loop begins. For 
example, if the independent variable is specified as RANGE, this range will 
be either total range at the end of the iteration loop or range to be flown 
in the initial segment of the iteration loop depending on the input option 
used in the initial segment. In any case the independent variable must be 
one of the parameters input in the segment data for the initial segment in 
the iteration loop. 

Data for Option 22. , 

The following data statements are used by Option 22; information in italics 
is optional. Note that the numbers a and 3 are decimal numbers and must 
contain a decimal point. The numbers A and B are integers. 


START loop at heginnirig of SEGMENT (A) 


FINISH loop at end of SEGMENT (B) 


NUMBER of iterations = a, 

DESIRED value of dependent var^hte - 3> 
END 

Example Input for Option 22 . 


WEIGHT I 

RANGE > is to he vari.ed 
MACH 1 


WEIGHT 

RANGE 

MACH 


is to he satisfied 


/ SEGMENT (27), OPTION (22) 

/ START LOOP AT SEGMENT (4) WEIGHT IS TO BE VARIED 

/ FINISH LOOP AT END OF SEGMENT(26)— — R/U^GE IS TO BE SATISFIED 
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/ NUMBER OF ITERATIONS = 3. 


/ DESIRED RANGE = 1500., NAUTICAL MILES 
/ END 


The default value for the number of iterations is 3. Note that if specified 
the number of iterations must be a decimal number. 

In the above example Segment (4) could, for example, utilize 0ption(13) to 
define a neiv value of weight and Segment (26) could utilize Option (6) to 
cruise to a specified weight. The initial value of weight in Segment (4) 
will be varied by Option (22) until the total range at the end of Segment (26) 
is 1500 ± 1.5 nautical miles or the number of iterations exceeds 3, whichever 
occurs first. 
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10. PRINT CONTROL INDICATORS 


IPDICN 

IPSEGi 

IPATMD 

IPTZIN 

IPENDS 

IPENDR 

IPTZXX 

IPRGSV 

IPDGPC 

AOUT 


Input dictionary print control; 

= 0, omit print 

= 1, print input dictionary in alphabetic and niimeric 
order 

Array indicating the final printout segment order and 
selection 

Atmospheric tabular print indicator: 

= 0, omit print 

= 1, print tabulation of atmospheric properties as a 
function of altitude 

Controls initial TERM array print 
= 0, omit print 

= 1, print term array from FILLTZ as each mission 
segment is processed 

Scaled propulsion system data print control 
=* 0, omit print 

= 1, print scaled propulsion data 

Raw (input) propulsion system data print control 
- 0, omit print 

= 1, print raw propulsion data 

Controls TERM array print as each segment is confuted: 

= 0, omit print 
= 1, print TZi 

Controls RG array print as each segment is computed 
= 0, omit print 
= 1, print RGj_ 

Generalized performance print control 
= 0, omit detailed print 
= 1, provide detailed print 

A five digit compressed term to control the printers; 

1st digit controls CHECK data output 

2nd digit controls PATH data output 

3rd digit controls CLIMB data output 

4th digit controls XCEL data output 

5th digit controls RUN and GOGO data output 

If AOUT is zero, the printing control will take on the value 
of the previous mission segment. The digit ”1" indicates 
DO NOT print out the data computed in the corresponding 
Subroutine; *'2" indicates print out all data computed in the 
corresponding subroutine 
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11. ALPHABETIC LIST OF INPUT AND INTERNAL VARIABLES 
11.1 A’s 


ACG 


Internal variable computed in CHECK 


ACV 


Internal variable computed in CHECK 


AERODYNAMIC MODE Optional way o£ speci^ing aerodynamic case number 
(CASE) 

AG Ground roll body angle 

AINLET Engine inlet area (not required if ramjet is to be scaled) (ft^) 

ALMXTO Angle of attack for take-off C^ MAX if internal table 

not to be used 

AM Initial Mach number for con 5 )utation of flight envelope 

AMACH Mach number array used in ML/D calculation 

A>IAXG Maximum ground rotation angle 

AMAXTJ Maximum Mach number for turbojet operation 

AME End Mach number for Mach-altitude paths on performance maps 

AMF Final Mach number for performance maps . Also end of climb 

Mach number FLITEN climb path history. If AMF=0.0, end 
of climb Mach number is determined from maximum range factor 

AMI Initial Mach number for performance maps 

AMINRJ Minimum Mach number for ramjet operation 

AMLIM Maximum allowable Mach number in FLITEN searches 


AMSWl First Mach number for propulsion system switch when CASEN=4. 

AMSW2 Second Mach number for propulsion system switch when 

CASEN=4. 


AMS 

AMTO 


Start Mach number for Mach-altitude paths on performance 
maps 

Initial Mach number for take-off 


AN Initial variable set in CHECK 

AOUT A five digit compressed teim to control the printers: 

1st digit controls CHECK data output 
2nd digit controls PATH data output 
3rd digit controls CLIMB data output 
4th digit controls XCEL data output 
5th digit controls RUN and GOGO data output 


If AOUT is zero, the printing control will take on the value 
of the previous mission segment. The digit "I" indicates 
DO NOT print out the data computed in the corresponding 
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AOUT Ccontd) subroutine; ”2" indicates print out all data cojnputed in 
the corresponding subroutine. 

AREFF Not used. 


ATABXX 

AIMOS 

AVOID 

AVOIDl. and 
AVOID2 


AVSFC 


Aerodynamic table XX 

Internal program array containing atmospheric properties 
at given altitude 

Compressed term 10 on First Mission Card 

Two digit terras to be used only when mission data change 
cards are used. (These terms must be zero the first time 
the mission is computed). The first variable, AVOIDl, 
represents a segment number, ni, and the last variable, 
AV0ID2, represents another segment number n2* AVOIDl 
equal to non-zero value indicates that the conputations in 
mission segments ni through nz will give the same results 
as the corresponding mission segments of the previous 
mission; therefore, segments nj through n2 will not be 
recomputed but will use the values computed in the previous 
mission . 

Internal variable computed in GOGO and RUN 
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11.2 B's 


BAREWT 

BFOBEX 


Vehicle bare weight 

Flap span to exposed wing span ratio. Used in take-off 
and landing calculations. 



11*3 C’s 


CASE 

CASEN 

CCPXX 

CD 

CDDELT 

CDO 

CDPCHT 

CDPLG 

CDSTOR 

CDTANK 

CFOCAV 

CL 

CLATOT 

CLC 

CLCD 

CLI.MBQ 


Aerodynamic case indicator for selecting systems 1, 2, and 3 
Propulsive case indicator for selecting systems 1, 2, and 3 


Compacted integer used in UPACK. Consists of 5 characters: 
1st character = propulsion option 
2nd character = aerodynamic option 
3rd character = throttle setting 

4th and 5th characters = number of integration steps 
NOTE : CCPXX is the 6th word in the Tt array printed after 

each segment data block is processed. 


Internally computed drag coefficient 

Drag coefficient increment used in FLIMLD 

Drag coefficient at zero lift. Landing and take-off 
analysis 

Flat plate drag coefficient of parachute, landing and take-off 
analysis 

Flat plate drag coefficient, landing and take-off analysis 

Internally computed stores drag coefficient 

Internally computed tanks drag coefficient 

Average flap chord to wing chord ratio, landing and take-off 
analysis 

Internally computed lift coefficient 

Linear lift curve slope per degree, landing and take-off 
analysis 

Internally computed lift coefficient 

Internally computed lift/drag coefficient 

Dynamic pressure for constant climb. If set to zero, 
dynamic pressure of End Point (AME, HE) will be used. 
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CLMAX 

CLMTAB 

CLMXTO 

CLO 

CLZZ 


Maximum cruise lift coefficient when INDCU4«0. (When 
INDCUl - 1. is found from CLMTAB in the form 

Table of as function of (M, CASE) when INDCLM * 1 

Maximum C. for take-off computations if internal tables 
are not to be used 

Lift coefficient at zero angle of attack. Landing and 
take-off analysis 

Internally computed maximum lift coefficient in FLIENV 



11.4 0»s 


DAC 

DATE 

DELAM 

DEiH 

DEurr 

DFK 

DMC 

Ditr 


Internally eorrputed aoeeleration in XCEL 
Date of a particular calculation 

Mach nuiri>er increaent in generalized perforwnce calcu* 
lations 

Increnental altitude in generalized perforaance calculations 

Weight decreaent or increaent in generalized parfoiaanca 
calculations 

Not ueed* Neferenoed in LANDIN 
internally eonputed Maoh number; need in PATS 
JHtemally eoiiputed drag foroe 
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END 


11.5 E's 


Signifies end of a data block 


ENDPTH 


End of path indicator for E climb; ' 

T , -rl , . Don/ 1 exceed terminal, a.ltitude t " 

= 0, terminal arc along final constant ^energy line 

» 1, Don't exceed terminal Mach number 


■ M “/"'V’’ / ii .1 

ETAB03. . . 

Engine data tables used in dual engine propulsion option 
: Secti-pn 6,v2.;3i for descra.pt ion ;, ^ .- f, - i i;,; 

ERRFP 

Undefined error indicator in PATH, internal 

ERRH 

Altitude error in PATH, internal 

:y.- .• ,■ -.v:* . f= 'i. ' 

ERRM 

Mach number error in PATH, internal 

.. y -.'r ■■ .i . 

Internal thrust error indicator, ENDAT and CMAX dumped 
from CRCOMP 

ElO 

Ell 

Internal error indicator, dunped from CRCOMP 


^ 7J-* 
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11.6 F*s 


FF 

FFLMIN 

FFR 

FFRMIN 

FLAPDC 

FLAPDT 

FN 

FNO 

FNMIN 

FNL 

FPOL 

FUEL 


Reserve fuel fraction (.05 for Mil C-5011A rules) 
Minimum fuel flow (un factored) per engine 
Intemdlty computed total fuel floiD rate 
Internally aonputed total fuel flou) rate JnCntmum 
Flap deflection, landing 
Flap deflection, take-off 
Internally computed net thrust per engine 
Internally set thrust in FILLPM and FLITOL 
Internally computed minimum total thrust 
Landing thrust ^ internally set 

Fraction of initial payload dropped or off-loaded during 
a mission 


FUELl 

FUEL2 

FUELS 

FUELIN 

FUELTl 

FUELT2 

FUELT3 

FUELTO 

FUELTK 

FUELUS 

FWC 


First fuel value, buddy refuel 

Second fuel value, buddy refue 

Third fuel value, buddy refuel 

Internally computed total internal fuel 

Internally computed fuel weighty tanks I 

Internally computed fuel weight , tanks 2 

Internally computed fuel weighty tanks 3 

Total fuel, including tank fuel at take-off 

Internally computed total tank fuel 

Internally computed total usable fuel at take-off 

Internally computed fuel weight change in CLIMB 
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11.7 G»s 


GLEVEL 

GLIMIT 

GRNDFL 

GRNDFT 


Number of 'g's' pulled in perforaence contour Mp 

calculations 

Limiting *g's' in radius of turn and tine to turn 
contour map 

Ground roll friction coefficient, landing 
Ground roll friction coefficient, take-off 
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HE 

HF 

HG 

HI 

HO 

HOVSFC 

HOVT 

HOVTMl 


End altitude on Mach> altitude path 
Final’ altitude on contour maps. 

Specified al^tudes in arbitrary atptosphere 

Initial altitude oh contour m^s 

Initial altitude on Mach -altitude path 

Not useds woe in OPT 2 7 

Not used; was in 0PT17 

Not usedj Uae in OPT27 



11.9 '^ V s 



IA40X« IA41X Number of first independent. Veriebl^ yeluef three- 
IA48X dimensional aerodynamic^^ ■ j 

IA40Y, IA41Y Number o/f S'e^d lnde^ vaiiiei in three- 

IA48Y dimensional aerodynamic tables 

.J ■ l'/ ■' ‘ I H ? 

ICLAX Table size imdicat^r, fox si .. . . .j V,- 

IClay ^ ' Lix ‘ ^ 


ICLZX Table size indicator for Ci 

ICLZY . Q . ,v, 

ICLMX Number of Mach number values., in CLMTitf . 

ICLMY Number of CASE values in CLMTAB 

IDRSPl Drop store pylon indicators: 

ISRSP2 * 0, Don't drop 

IDRSP3 * 1 * drop 

lORTPl Drop tank pylon indicators: 

IDRTP2 * 0, don't drop 

IDRTP3 ■ 1» drop 


lEOP 

IE03X 

IE04X 

IE05X 

IE06X 

IE06X 

IE07X 

IE08X 

IE09X 

lEOlOX 

INDE03 

INDE04 

INDE05 

INDE06 

INDE07 

INDE08 

INDE09 

INDEIO 

INDCLA 

INDCLZ 


Table size indicators for combined engine simulation data 


Table usage indicators for combined engine data 


Table usage indicator for Cl„ 
Table usage indicator for 
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INDAOl... Aerodynamic table indicators 

INDA12, 

INDA40, INDA41 
INDA60 

INDAER Clean vehicle aerodynamic indicator 

- 8, general model 
=9, derivatives 


INDATM 


INDCLM 


INDGPC 


Atmosphere indicator 

= 3, use original NSEG model 

- 9, use arbitrary model atmosphere 

Controls CL^ax selection: 

= 0, CLjnax given by CLMAX value 
~ 1> CLjnax ” CASE) : CLhfTAB 

Flight envelope calculation control indicator 


INDMAP 


INDMAP 

(Contd.) 


INDMHP 


INDTOA 

INDTOl 
. . . . INDT06 


Selects which performance contour map is to be drawn in 
current data block: 

= 0, no maps 
= 1 , E map 
= 2 , E/W map 
= 3> L/D map 
=4, Range factor map 
= S, thrust map 
^ 6 , drag map 

= 7, §FC map 
8, W map 
= 9, energy map 
10, L/(T-D) map 
= 11, turn radius map 
= 12, time to turn map 

Controls contours to which tangency criteria is applied 
(see IPATH) : 

= 1, use max [E] condition 
= 2, use max. [E/W] condition 
= 3, use max. [L/D] condition 
=4, use max. [range factor] condition 
* 5, use max. [L/(T-D]] condition 

Nonzero specifies that CLMXTO and ALMXTO will be used for 
take-off and landing computations in lieu of internal tables 

Propulsion table indicators 


INDTRN Turn calculation Cl selection indicator: 

= 0, turn at constant velocity 
= 1, turn using Cj^^ 

INDZLT Selects contour levels on performance maps : 

= 0, program computes levels 
- 1, levels specified in ZLTAB 
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IPATH 


IPATMD 

IPDGPC 
IPDICN 
IPENDR 
I PENDS 

IPLFl 

IPLF2 

IPLF3 

IPRGSV 

IPSEG 

IPXX 


IPTZIN 

IPTZXX 


Controls tangency condition for Mach*altitude path inscribed 
on performance contour maps: 

= -1, no paths 
= 0, ail paths 

= .1, energy tangency climb 
= 2, altitude tangency climb 

= 3, Mach tangency climb 

= 4, linear M-h climb 

= S, constant q climb 

Generate and print atmospheric table indicator* Controls call 
to Overlay 3.1 in MISI(X*4. 

Detail print indicator for flight envelope calculations 
Data base directory print control 
Print raw propulsion system data indicator 
Print scaled propulsion system data indicator 


Indicator to control allocation of tank fuel to payload: 
= 0, don't include in payload 
= 1, include in payload 

Print the RG array from each mission segment indicator 


Array giving mission segments to be printed and their 
order at printing 


Array of 
IPXX(l) , 
( 2 ) 

(3) 

(4) , 

(5) , 

( 6 ) , 

(7) , 

( 8 ) , 

(9) 

( 10 ) 


print control indicators: 

controls RERAN, CRCOMP, and CHECK print if .TRUE. 


controls PATH print of .TRUE, 
controls CLIMB print if .TRUE, 
controls XCEL print if .TRUE, 
controls RUN print if .TRUE, 
controls GOGO print if .TRUE. 


NOTE THAT in FLITEN calculations the IPXX array is set 
internally to .TRUE. 


Print term array as each mission segment data is processed 
indicator 


Print term array as each mission segment is started 
indicator 
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ISEG 

ISEC 

ITOIX,. .. 

. . . IT06X 

ITOIY, 

. . . IT06Y 

IT02W 

IT04W 

ITZ 


Internal variable uhioh oounte eegment nwrber, (ElfSEG) 

Internally bdniputed ^aeOond ixalue in ealle tt> ’CLOCK 

Number of values for fitst va:ri^bl'e in three^ or' -four 
dimensional propulsion tables 

Number of values for second variable in three or four 
dimensidhil pro^ulisioh tables ' " , > 

Number of values for third 'variable ih four dimensional 
propulsion tables 

Internally oonputed location for yoeition of segment in 
TZ array 



JDRPSl 


ii.io j*s 

NuniJBT of 8tora pairs type 1 to be dropped 
JDRPS2 Itunter of store pairs type 2 to be -dropped. , 

JDRPS3 fharber of store poire type 3 to be dropped 

JDRPTl Nwiber of tank pairs type 1 to .be <3^^ 

JDRPT2 Umber of tonk pairs type 2 to be cirppp^d 

■'•-■i. . '.0. »«.5 ■ ;V 

JDRPT3 Nwtber of tank pairs type 3 to be dropped 
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11.11 L*s 


LOCLM 

LOCNG 

LOCRG 

LOCSAC 

LOCTFF 

LOCTZ 

LTTIOX 

LTTIOY 

LTTIOW 

LTTIOZ 

LTTllX 

LTTllY 

LITllW 

LTTllZ 


Internal^ location of CLMTAB in TAB array 

Internal, location of RG array in C array 

Internal, etarting location of aerodynamic tdblee in 
TAB array 

Internal, etarting location of propulsion tables in 
TAB array 

Internal, location of TZ array in C array 
Not used 
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11.12 M*s 


MAXSEG 

MDROPl 

MDROP2 

MDROP3 

MIN 

NTTABOl 

MTAB02 

MTAB03 


Maximum number of mission segments 

Nwtber of tank poire type 1 dropped^ 
Nuniser of tank pairs type 2 dropped, 
Number of tank poire type 3 dropped, 

Intematty set minutes from CLOCK 

Mission segment data tables 


intemat 

internal 

internal 
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11.13 N's 


NAM 

NATIO 

NATH 

NCASE 

NDROPl 

NDROP2 

NDR0P3 

NH 

NHG 

NPAGE 

NSDROP 

NSTORl 

NSTOR2 

NSTOR3 

NTANKl 

NTANK2 

NTANK3 

NUMCON 

NUMENG 

NUMSEG 

NWT 


Number of Mach numbers in generalized performance calculations 


Number of stores type 1 dropped, internal 

Number of stores type 2 dropped, internal 

Number of stores type 3 dropped, internal 

Number of altitudes in generalized performance calculations 

Number of altitudes employed in arbitrary atmosphere model 

Page aountery internal 

Number of store pairs type 1 
Number of store pairs type 2 
Number of store pairs type 3 
Number of tank pairs type 1 
Number of tank pairs type 2 
Number of tank pairs type 3 
Number of contours in performance map 
Number of engines 

Internally computed mieeion segment numbers set to 1 in 
INISEG 

Number of take-off weights for climb path 
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Zntmxnatty 99t mieaian cpticn indCoator 


OPTION 



11.15 P*s 


PAIRS 

PAYWT 

PLWINT 

PM 

PQl 

PQ2 

PQ3 

PQ4 

PQ5 

PQ6 

PQ7 

PQ8 

PQ9 

PQIO 

PS 

PSET 


Number of tank pairs carried in simple tank option. 

See ENERGY, 0PTXX3 

Internally oomputedj total payload 

Internal payload weight 

Array of specified pressures at base of each layer in 
arbitrary atmosphere model; in pounds per square foot. 
If pressures are free, omit. 

Internal t computed in OUTXXX, SEGOUT 

Internal^ computed in OUTXXX^ SEGOUT 

Internal^ computed in SEGOUT 

Internal^ computed in SEGOUT 

Not used 

Not used 

Not used 

Not used 

Not used 

Not used 

Power setting, propulsion 
Throttling option: 

= 1, maximum available dry, wet or maximuiB power 
=2, minimum available 
s 3, at specified power setting 
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11.16 Q's 


QLIM 

QS 


Maximum allowable dynamic pressure 

Internal t ■product of dynamic pressure and wing areas 
computed in CHECK j PATH;, CRCOMPy EDMHFTj etc. 



11.17 R»s 


RANGE 

RANGEl 

RANGE2 

RANGES 

RC 

RCM 

RCMAX 

RCSC 

RETWGT 

RFAC 

RN 

RXC 


Internally computed range ^ EN^G 

Three, range values for buddy refuel mission segment option 

Internal^ robe of climb; computed in FAIH 
Internal range increments FLIEEV 
Internals used in OFTXXZ 
Rate of clinib at service ceiling 
Internals return weight 

Internals Breguet range factor; computed in CBECKs CECOMP 
Internals eurnned range FLIMHC 

Internals range covered in time BT; computed in CLIMB 
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11.18 S’s 


SFC 

SFCFAC 

SPCHUT 

SPTYPl 

SPTYP2 

SPTYP3 

SPLG 

SQ 

SREF 

STWRET 

STWDRP 

STYPl 

STYP2 

-STYP3 


iHterridt, '6peai.fi^o fiiet adnsimption^ ENDAT 
Factor applied to specific fuel consumption 
Flat plate area of drag chute 

Internal i number of store pylons type 1 left^ SACSP^ STORES 
Internal i number of store pylons type 2 left^ SACSPj STORES 
Intemalj number of store pylons type 3 left, SACSP, STORES 

Flat plate area of landing gear 


Wing reference area 


Internal, number of stores type 1 left, SACSP, S. ' 

Internal, number of stores type 2 left, SACSP, Sl .a. 

Intemal, number of stores type 3 left, SACSP, STOR^f^ 
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TAG 

TCMAX 

TERM 

TEWT 

TFWT 

TIME 

TITLE 

TM 

TMINRJ 

TMINTJ 

TOL 

TOTSWD 

TOTTED 

TOTTWD 

TPWT 

TPTYPl 

TPTYP2 

TPTYP3 

TRATIO 

TSIZEF 

TSWT 

TSWDRP 

TSWDRP 

TSWRET 


11.19 T's 

Internal^ time to acoeie-pate , XCEL 

Internal j time variable OPTXX2y OPTXX4 

Internal^ array of partiaular miseion segment data 

Weight of empty fuel tank pair, simple tank option (old NSEG) 

Weight of fuel in tank pair, simple tank option (old NSEG) 

Internal j mission timoy ENSEG 

Title card, read in INSEG 

Array of temperatures used in arbitrary atmosphere 

Minimum ramjet throttle setting for combined engine option 
Minimum turbojet throttle setting for combined engine option 
Internal y tolerance set in CLIMB 

Internal y total store and pylon weight droppedy STORES 
Intemaly total tank fuel dropped or aonsumedy TANKS 
Intemaly total weight of tank and pylons droppedy TANKS 

Intemaly number of tank pylons type 1 lefty SACTPy TANKS 

Intemaly number of tank pylons type 2 lefty SACTPy TANKS 

Intemaly nurrber of tank pylons type S lefty SACTPy TANKS 

Wing taper ratio, landing and take-off aerodynamics 

Engine sizing factor 

Intemaly total tank and store system weighty TNKSWT 
Internal, total store system weight dropped 
Intemaly total store system weight returned 
Intemaly total store system weight returned 
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TTABXX 


TTCC 

TTR 

TTYPl 

rryp2 

TTYP3 

TZ 

TZ13 
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Propulsion table number XX 
Internal variable in CLIMB 

Internal^ time to aruiaes RUN; uaed in several options 

Internal^ number of tanka type 1 leftySACTP^ TANKS 

Internal^ nunber of tanks type 2 left^SACTP, TANKS 

Internals nunber of tanks type 3 leftjSACTPj TANKS 

Internals array of mission segment data* Contains all 
segments and is aggregate of TERM arrays 


I 


11,20 V’s 


VN 

vx 

V50 


Internals velocity knotSj CHECK 

Internals velocity in feet per eeconds CHECK 

Internals velocity at SO foot obstacles take-off analysis 
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11.21 w*,s . 


w 

Initial weight 


W1 

Intemdl, weight FLIMSC 


W2 

Intexnta%:t weight. FhlMHC 


W50 

IntBxnat, weight at, SOr foot obetaole, take-off cpuityBia 

WFF 

Final weighty endurance and range performance 

WGTO 

Initial weight 


WGTFTl 

Weight of fuel in tank pair type 1 


WGTFT2 

Weight of fuel in tank pair type 2 


WGTFT3 

Weight of fuel in tank pair type 3 


WGTSl 

Weight of store pair type 1 


WGTS2 

Weight of store pair type 2 


WGTS3 

Weight of store pair type 3 


WGTSPl 

Weight of store pylon pair type 1 


WGTSP2 

Weight of store pylon pair type 2 


WGTSP3 

Weight of store pylon pair type 3 


WGTTl 

Enpty weight of tank pair type 1 


WGTT2 

Empty weight of tank pair type 2 


WGTT3 

Empty weight of tank pair type 3 


WGTTPl 

Weight of tank pylon pair type 1 


WGTTP2 

Weight of tank pylon pair type 2 


WGTTP3 

Weight of tank pylon pair type 3 


WGTUFl 

Weight of usable fuel in tank pair type 

1 

WGTUF2 

Weight of usable fuel in tank pair type 

2 

WGTUF3 

Weight of usable fuel . in tank pair type 

3 
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WI Initial weight, endurance and range performance 

WINCID Wing incidence used in take-off and landing analysis 

WL Weight at start of landing; also control landing analysis 

in FILLPM 

WO Initial take-off weight; also contains take-off analysis 

in FILLPM 

WSLEFC Wing leading edge sweep angle; used in landing ant take-off 

analysis 

WT Weight used in general performance calculations 

WTKEP Internal set to initial segment weight in 0PTXX7 and other 

options 

WUTO Warm up and take-off fuel weight allowance 

WUTOFU 
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11.22 X's 


XAl 

XA2 

XA3 

Internal variables in FLYMHP, FLIXCP^ OPTXX2 


XAA 

Internal 3 controls INDMHP setting in CMAX» Control maximum 
R/C or dh/dW in OPTX14 

XAO 

Internal dM/dh in FLYMHF^ FLIXCP, OPTXXl 


XBO 

XBl 

Internal^ used in EMAX^ FLIENV^ FLIXEP, FLIXCP, XBO 
XBl AMLIM in OPTX24 

QLIM, 

XFLEN 

Internal indicator. Controls call to FLITEN in MISION, 
Set in ENSEG 

XISO 

Specific impulse at take-off 


XiMAXG 

Maximum ground roll 


XMLOD 

Internal^ set in RNSEG 


XMSEG 

Internal i used in INSEG^ set in RNSEG 


XNZZ 

Internal, error indicator set in BERAN, CHECK, CLIMB, 
OPTXXl, OUT XXX, etc. ' 

CRCOMP, 

XQ12 

Internal, used in 0FTX14 


XQAl 

Internal, altitude variable GOGO OUTXXX 


XQAE 

Internal, altitude variable GOGO, OUTXXX 


XQBD 

Internal, indicator controlling CHECK Option and used 
EMAX, FLIMHC, FLIXCP, GOGO 

in 

XQFP 

Internal constant, FLIMHP, FLIMHC, FLIXEP, FLIXRP, FLIXCP, 
OPTXXl, OUTXXX, etc. 

XQRl 

Internal, range variable GOGO, RUN 


XQR2 

Internal, used in RUN 


XQR3 

Internal, used in RUN 


XQSW 

Aerodynamic reference area 
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XQTl 

XQT2 

XQT3 

XQWl 

XQW2 

XQWE 

XSTEP 

XTT 


Internal, time variable GOGO, RUS 
Internal, used in RUN 
Internal, used in RUN 

Internal, weight variable in GOGO, OUTXXX, RUN 
Internal, weight variable in GOGO, OUTXXX, RUN 
Internal, minimum weight in CLIMB 

Internal, limiting altitude ahange in FLXMBP* Aleo used in 
FLIMBC, FLIXEP, FLIXRP, FLIXCP, GOGO, INIMHP, OPTXXl, eta. 



11.23 Y's 


YZA 


IntemaZj vooctabZe used in CHECKS CLIMBj PATHj eta. 
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11.24 Z's 


ZLTAB Table of "Z" levels, contour plotting, used when INDZLT « 1. 

ZMAXIO Minimum L/(T-D); map 10 

ZMAXll Maximum "Z" contour, map 11 (turn radius) 

ZMAX12 Maximum "Z'* contour, map 12 (time to turn) 

ZXS Internal^ nurribev of tanke dropped^ aimpte option, ENERGY 

ZX6 Internal, used in PATH, RUN, ZX6 - 0.0 in SEGINI 

ZX7 

ZXS Internal, set to WGTO in SEGINI 

ZX9 Not used 

ZXIO Internal, weight term in ENERGY, Set to total non^fuel 

weight dropped in OPTX12. Set to PAYWT in SEGINI 

ZYl 

ZY2 

ZY5 

ZY4 

ZYS Not used 

ZY6 

Til 

Ti^ 

Ti9 

ZYIO 
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12. FLIGHT PATH INEQUALITIES 


AMAXTJ 

AMINRJ 

AMLIM 

CLMAX 

CLNfTAB 

FNMIN 

FFRMIN 

GLEVEL 

QLIM 

TMINRJ 

TMINTJ 


Maximum Mach number for turbojet operation 
Minimum Mach number for ramjet operation 

Limiting flight Mach number 

Limiting Cl value when INDCLM - 0 

Table of Cr as function of (MACH, CASE) when 
INDCLM = 1^^^ 

Minimum thrust 

Minimum fuel flow 

Load factor 

Dynamic pressure limit 

Minimum ramjet throttle setting for combined engine option 
Minimum turbojet throttle setting for combined engine option 



13. GENERALIZED PERFORMANCE ANALYSIS 
Two types of generalized performance calculations are available in NSEG; 

1. Contour maps of specified performance function 
with or without superimposed climbs, descents, 
or accelerations 

2. Organized sets of mission segments for a range 
of weights , Mach numbers , and altitudes 

Generalized performance calculations are controlled from subroutine FLITEN 
which is called directly from the Main program of Overlay (3.0) MISION when 
the internal variable XFLEN is .TRUE. Contour maps are prepared in FLIMAP 
Overlay (3.4), and organized missions are prepared in FLIRAN Overlay (3.3), 
Figure 17 presents a flow chart for subroutine FLITEN. 

Each map generated requires a self-contained block of data which terminates 
with an END card. Vehicle data and contour map data, such as map limits, may 
be altered in any map data block. All map data blocks are contained in the 
GENPC data region. Figure 17. 


13.1 Contour Plotting Data and Options 

A set of point calculations (i.e., vehicle capability at given flight 
conditions) are carried out over a two-dimensional array of Mach/altitudes, 
Mi, hj . The resulting matrix of capabilities, F^j is then supplied auto- 
matically to the CONPLOT routines, and the contours of the function F^ in 
the Mach/altitude plane are obtained in the form of CALCOMP, Houston plotter, 
or CRT display device output. At the present time twelve functions, F^ to 
f 12, may be output in contour form. Each contour plot and input is described 
below. 

AME End climb Mach number 

AMF Highest Mach number on contour grid 

AMI Lowest Mach number on contour grid 

AMS Start climb Mach number 

ENDPTH Terminal maneuver indicator: 

=-l, don’t exceed terminal altitude 
= 1, allowed to exceed terminal altitude 
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GLEVEL Number of 'g* pulled by vehicle 

he End climb altitude 

HF Highest altitude on contour grid 

HI Lowest altitude on contour grid 

HO Initial climb altitude 

INDGPC Generalized performance calculation indicator. Must be set 
to INDGPC « 8 for maps 

INDMAP Map indicator: 

■ 0, no map 
s i, produce map i 

INDZLT Contour level indicator: 

« 0, Internally computed by NSEG 
a 1, specified directly by user 

IPATH Superimposed path indicator: 

a-1, no paths 
a 0, all paths 
a 1 , I climb 

a 2, h climb 

a M climb 

a 4, linear M-h climb 
a 5f constant dynamic pressure climb 

IPDGPC Controls print out of grid details: 
a 0, don't print 
a 1, print 

NAM Number of Mach numbers in grid, (Max. = 21) 

NH Number of altitudes in grid (Max. » 21) 

NUMCQN Nuober of contour levels to be used 

ZLTAB Table of contour levels to be defined on map when INDZLTal 
In addition, the appropriate aerodynamic and propulsive data must be specified. 
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13.1.1 


Specific Energy Time Derivative, E (INDMAP=1) 

The specific energy time derivative is computed according to the expression 

E(M,h) = CT-D)V/W 


where 

E - energy total time derivative 

T “ thrust obtained at a specified power setting or at T=D. Wet, dry, 
or maximum power options are available 
D - drag computed for a specified load factor 
V - flight velocity 
W - aircraft weight 


13.1.2 Specific Energy/Fuel Flow, E/m, (INDMAP=2^ 

The E/ra contour presents the specific energy time derivative over the fuel 
f 1 ow map . Since 

• /• _ dE/dt _ ^ 
dm/dt dm 

The point calculation employed is 

E/m = (T-D)V/(Wra) 
where m is the fuel flow rate. 


13.1.3 Lift/Drag, L/D, (INDMAP=3) 

Lift/drag contours present a measure of the airplane’s aerodynamic efficiency. 
The L/D maps indicate its range capability in unpowered flight and partially 
reflect the cruise range capability. Maps can be produced for any specified 
load factor. 


13.1.4 Range Factor, Rp, (INDMAP=4) 


Range factor contours present a measure of vehicle cruise range capability. 
Maps are produced for level flight with thrust equal to drag at a specified 
aircraft weight. 


R = r JL. f L-) 
F '^SFC-^ 


where SFC is the specific fuel consumption. 


13.1.5 Thrust (INDMAP=5) 

The thrust map is available as a device for examining thrust input data 
or the thrust component of other mapped functions. The map can be obtained 
for wet, dry, maximum, or throttled power setting. 


120 



13.1.6 Drag (INDMAP=6) 


The drag map provides a device for inspecting drag data input or the drag 
component of any other map. Drag maps are produced for a specified load 
factor. 


13.1.7 Specific Fuel Consumption, SFC, (INDMAP=7) 

Specific fuel consun^tion maps are provided as a data input inspection device 
or as an aid to visualizing the specific fuel consumption component of other 
maps. Maps may be obtained for wet, dry, maximum, or throttled power settings. 


13.1.8 Fuel Flow Rate, ra, (INDMAP = 8) 

The fuel flow maps are provided as a data input inspection device or as an 
aid to visualizing the fuel flow component of other maps. Maps may be 
obtained for wet, dry, maximum, or throttled power settings. 


13.1.9 Specific Energy (INDMAP=9) 


The specific energy map 

E = h + V^/2g 

is provided as a user's convenience in visualizing the trajectory points 
between constant energy lines and any other set of contours. 


13.1.10 Lift/ (Thrust- Drag) , L/ (T-D) , (INDMAP=10) 

The lift/ (thrust-drag) contours are useful for determination of maximum reinge 
powered flight. Assuming that maximum range flight occurs at small flight 
path angles with L = W 

^ / T^ 

Therefore, the energy- like approximation to maximum range flight occurs when 
L/(T-D) is a maximum at each energy level. It should be noted that when 
T - D = 0, no energy gain is possible; therefore, this singular condition must 
be avoided. In NSEG the per cent excess of thrust over drag which is 
acceptable is the program input ZMAXIO. 
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13.1.11 Turn Radius (INDMAP = 11) 

Turn radius maps give a gross indication of aircraft's combat capability. 
Turn radius is computed by equating the aircraft's lift capability in 
steady state or decelerating flight using the e:q)ression 



where is determined so that 

1. INDTRN = 0 

2. thrust equals drag for steady state flight 

3. INDTRN == 1 


4. 


Cl = 


Ct for minimum instantaneous turn radius 
%iax 


The maximum turn radius to be plotted is given by ZMAXll. 


13.1.12 Time to Turn, (INDMAP=12) 

Time to turn through 180® is presented as a supplement to the turn radius 
When the minimum instantaneous turn radius calculation is employed, the maps 
do not give a true time to turn. They merely indicate how long a time the 
aircraft would take to turn if it could maintain its current turn rate* When 
steady state turns are considered, true time to turn is obtained which will 
frequently be rruch longer than is required for a decelerating turn. 

The maximum time to turn 180® plotted is given by ZMAX12. 


13.2 Organized Mission Segments and Flight Envelope 

These calculations are retained in the original NSEG program form. General 
data required is 

AMLIM Maximum Mach number 

CLMAX Maximum lift coefficient 

QLIM Dynamic pressure limit 

Appropriate propulsion and aerodynamic data must also be specified. 
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13.2.1 Flight Envelope Calculation (IIJDGPC = 1) 

AM Initial Mach number for computation of flight envelope 

DELAM Mach increment for flight envelope 

DELWT Weight increment for computation of several flight envelopes 
NWT Number of weights for confutation flight envelope 

WT Initial weight for computation of flight envelope 

13.2.2 Climb Path History, (INDGPC=2) 

AMF ' End of climb Mach number (if zero, RFMAX will be used) 

DELWT Incremental weight for climb performance 

NWT Number of weights to compute climb path 

WT Initial take-off weight for climb performance 

WUTO Warm-up and take-off fuel 

13.2.3 Endurance vs. Weight at Several Altitudes, (INDGPC = 3) 

DELH Incremental altitude 

H Initial altitude for endurance computation 

WFF Final weight for endurance 

WI Initial weight for endurance computation 

The number of weights, NWT, is computed internally and is less than or equal 
to 9. 

13.2.4 Optimum Cruise Climb vs. Weight at Various Mach Numbers (INDGPC=4) 

AM Initial Mach number for cruise climb 

DELAM Incremental Mach number 

NAM Number of Mach numbers 

WI Initial cruise weight 

WFF Final cruise weight 
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The number o£ weights, NWT, is internally computed and is less than or equal 
to 9. 

13.2.5 Constant Altitude, Constant Mach Cruise at Various Weights, (INDGPC = 5) 

AMI Initial Mach number for constant altitude cruise 

DELH Incremental altitude 

DELAM Incremental Mach number 

NH Number of altitudes 

NAM Number of Mach numbers 

WFF Final weight for cruise 

WI Initial weight for cruise 

The number of weights, NWT, is internally computed and is less than or equal 

to 9, 

13.2.6 Take-Off and Landing Calculations (INDGPC = 6) 

Data has already been presented in Section 8. 

13.2.7 ML/D Calculation at Given Altitude, Various Mach 
Numbers and Cj^ Values, (INDGPC = 7) 

AMACH Mach number value table 
H Given altitude 

NAM Number of Mach numbers, less than or equal to 9. 

Note that ML/D is computed for Cj^ in the range 0 to 1.6 at .05 Cj^ intervals. 
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OVERLAY 0.0 


OVERLAY 1.0 
INISEG 

• Initialize 


NSEG II 


OVERLAY 2.0 
INSEG 

• Read Input 



Print Atmosphere 


Missions 


FIGURE 1. NSEG 


♦ MAIN CONTROL 


OVERLY 

\Y 3.0 

MIS] 

[ON 


1 

1 

OVERLAY 3.3 


OVERLAY 3.4 

FLIRAN 


FLIMAP 

■ Performance 

» Contour Maps 


• Control 
Calculations 


STRUCTURE 




FIGURE 2. BASIC NSEG DATA OUTPUT 





(a) 

(b) 

1. 

TITLE 


1. 

TITLE 

2. 

TABLE SIZES 


4. 

GENERAL PERFORMANCE 

3. 

VEHICLE DATA 


6. 

RUN CARD 

4. 

GENERAL PERFORMANCE 



5. 

MISSION 



6. 

RUN CARD 

. 




(c) 


(d) 

1. 

TITLE 


1. 

TITLE 

5. 

MISSION 


3. 

VEHICLE DATA 

6. 

RUN CARD 


5. 

MISSION 



6. 

Run Card 


(e) 


(f) 

1. 

TITLE 


1. 

TITLE 

2. 

TABLE SIZES 


2. 

TABLE SIZES 

3. 

VEHICLE DATA 


3. 

VEHICLE DATA 

4. 

GENERAL PERFORMANCE 


5. 

MISSION 

6. 

RUN CARD 


6. 

RUN CARD 

FIGURE 3. SOME POSSIBLE ADDITIONAL NSEG DATA BLOCKS 
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START 



RETURN 


Provides automatic propulsion 
switching on Mach number 


FIGURE 4. FLOW CHART, ENCASE 






START 


I 



• Write Scaled TTAB02 

FIGURE 5. FLOW CHART ENDAT 






WET POWER 


C 


INDT03 ^ 

- ■ I 

IPENDR « IPENOS « O}- 


C 


fPAGHED I 


IPENDR = 0 


3 


[write: “MAX> thrust CASEN = 2. RAW DATA” | 


T 

JAhmach i 

i ASCALE I 

± 


Writes Out TTAB03 


Scales TTAB03 


^ IPENDS * 0 ') 

± 


WRITE: "MAX. THRUST CASEN=2 SCALED DATA" 


< 


I ^ 1 

I AHMACH I 


INDTU4 


> 


■ Writes Scaled TTAB03 
=0 


■*o 


C 


IPENDR » 0 


>— 


i t^AGHED I 

i 


I WRITE: "SPECIFIC FUEL CONSUMPTION CASEN^Z, RAW DATA" 


I 

I pHfiACH"! 

I bscaleI 

IPENDS » 0 y 


Writes out TTAB04 


Scales TTAB04 
Y 


<D 


I PAGHED I 


WRITE: "SPECIFIC FUEL CONSUMPTION CASEN=1 , 
SCALED DATA" 


1 


I PHMACH I 

• Write Scaled TTAB04 

FIGURE 5. FLOW CHART ENDAT (CONTINUED) 


*© 
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MAXIMUM POWER 


WRITE: "MAX. THRUST CASEN=3, RAW DATA" 


AHMACH 


ASCALE 


• Writes Out TTAB05 


Scales TTAB05 


WRITE: "MAX. THRUST CASEN=3, SCALED DATA" 


AHMACH 


Writes Scaled TTAB05 




PAGHED 


WRITE: "SPECIFIC FUEL CONSUMPTION CASEN = 3, RAW DATA" 


CHMACH 


CSCALE 


Writes Out TTAB06 
Scales TTAB06 


WRITE: "SPECIFIC FUEL CONSUMPTION CASEN » 3, SCALED DATA" 


Writes Scaled 
TTAB06 



FIGURE 5. FLOW CHART ENDAT (CONCLUDED) 
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y 

RETJRN 


FIGURE 6. FLOW CHART ENSET (CONTINUED) 
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FREESTREAM DIRECTION 





LIFT 



PARALLEL TO WING 
UNDERSURFACE 



= SPILLED AIRFLOW/ 

^2 " A^/TAN5j 

^2 = - Pi)A2 

^SPILL = ^^2 ^SPILL = ^ ^2 

FIGURE 8. SPILLAGE LIFT AND DRAG COMPUTATIONS 
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COMPUTE TOTAL FREESTREAM 
ENTHALPY 

HO= H+VO*VO/50073.2 
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COMPUTE: 

SPILLED AIRFLOW (WASRJ) 
SPILLAGE LIFT AND DRAG 
DUE TO RAMJET SPILLAGE 



COMPUTE: 

TOTAL SPILLAGE LIFT AND DRAG 
RAMJET FUEL FLOW RATE (WFRAM) 

RAMJET TOTAL THRUST (FTOTRJ) 

RAMJET NET THRUST IN FREESTREAM DIRECTION (FAVRAM) 
TOTAL FUEL FLOW RATE (WFTOT) 

TOTAL AIRFLOW RATE (WATOT) 

TOTAL NET THRUST (FAVTOT) 

RAMJET THROTTLE SETTING (TSETRJ) 



TOTAL NET THRUST 
EQUAL THRUST 


YES 




























START 


CLEAN 

VEHICLE 



RETURN 





START 











Cj3=Min Cf). : > 0.0 

i ^ ■ 


S~Si'*'^D2'*'^D3 


RETURN WRITE: "Cp.=i, CASE=j” 


FIGURE 13. FLOW CHART, SUBROUTINE SACSN2 
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RETURN 




165 


START 



RETURN 
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SEARCH FOR MACH 
ALTITUDE POINT 
USING 

BASIS FOR SELECTION OF MACH ALTITUDE POIOT 

MAX 

E 

NIAX 

E/M 

MAX 

5,/D 

MAX- 

Range 

Factor 

MAX 

h/(fn-drt) 

Constant 

Energy 

IPATH=1 

INDMHP=1 

INDMHP=2 

INDMHP=3 

INDMHP=U 

INDMHP=5 

Constant 

Altitude 

IPATH=2 

INDMHP=1 

INDMHP=2 

INDMHP=3 

INDMHP-U 

INDMHP=5 

Constant 
Mach Number 
IPATH=3 

INDMHP=1 

INDMHP=2 

INDMHP=3 

INDMHP=U 

INDMHP=5 

IPATH=U 

Linear Mach Altitude Path 
Fly from (HO, AMS) -> (HE AME) 

IPATH=5 

Fly Constant Q Path 


If IPATH = 0, all five paths 
will be generated using the 
function specified by INDMHP 
as the basis for selecting 
each Mach altitude point along 
the path. 


If IPATH < 0 , no path will be 
generated. 


FIGURE lb , PATH FINDER OPTIONS OF NSEG MISSION OPTION 1, : 







FIGURE 18, FLOW CHART FLITEN 
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